Mon. Not. R. Astron. Soc. 000, [T]-?? (2011) Printed 24 February 2013 (MN WI^ style file v2.2) 



Satellite abundances around bright isolated galaxies 

Wenting Wang^'^'^, Simon D.M. White^ 

^Key Laboratory for Research in Galaxies and Cosmology of Chinese Academy of Sciences, M ax- Panck- Institute Partner 

Group, Shanghai Astronomical Observatory, Nandan Road 80, Shanghai 200030, China 

^ Max Planck Institut fur Astrophysik, Karl-Schwarzschild-Str. 1, 85741 Garching b. Miinchen, Germany 

■^Graduate School of the Chinese Academy of Sciences, 19A, Yuquan Road, Beijing, China, 100080 



o 



O: 

u 

6 

^' 

> 

O 

o 
p 

rn 

O 



'i 



Submitted to MNRAS 



ABSTRACT 

We study satellite galaxy abundances by counting photometric galaxies from the 
Eighth Data Release of the Sloan Digital Sky Survey (SDSS/DR8) around isolated 
bright primary galaxies from SDSS/DR7. We present results as a function of the lu- 
minosity, stellar mass and colour of the satellites, and of the stellar mass and colour 
of the primaries. For massive primaries (logM^/MQ > 11.1) the luminosity and stel- 
lar mass functions of satellites with log M^/M© > 8 are similar in shape to those 
of field galaxies, but for lower mass primaries they are significantly steeper, even 
accounting for exclusion effects due to the isolation criteria. The steepening is par- 
ticularly marked for the stellar mass function. Satellite abundance increases strongly 
with primary stellar mass, approximately in proportion to expected dark halo mass. 
For log Mj,/M0 > 10.8, red primaries have more satellites than blue ones of the same 
stellar mass. The effect exceeds a factor of two at log A'/*/Mq ^ 11.2. Satellite galaxies 
are systematically redder than field galaxies of the same stellar mass, except around 
primaries with logM^/M© < 10.8 where their colours are similar or even bluer. Satel- 
lites are also systematically redder around more massive primaries. At fixed primary 
mass, they are redder around red primaries. We select similarly isolated galaxies from 
mock catalogues based on the galaxy formation simulations of Guo et al. and analyze 
them in parallel with the SDSS data. The simulation reproduces all the above trends 
qualitatively, except for the steepening of the satellite luminosity and stellar mass 
functions with decreasing primary mass. Model satellites, however, are systematically 
redder than in the SDSS, particularly at low mass and around low-mass primaries. 
Simulated haloes of a given mass have satellite abundances that are independent of 
central galaxy colour, but red centrals tend to have lower stellar masses, reflecting 
earlier quenching of star formation by feedback. This explains the correlation between 
satellite abundance and primary colour in the simulation. The correlation between 
satellite colour and primary colour arises because red centrals live in haloes which are 
more massive, older and more gas-rich, so that satellite quenching is more efficient. 

Key words: galaxies:abundanccs-galaxies:cvolution-galaxies:luminosity function, 
mass function-galaxics:statistics-cosmology:obscrvations-cosmology:dark matter 



1 INTRODUCTION 

Clustering studies provide insights into the formation and 
evolution of galaxies that complement those coming from the 
joint distribution of intrinsic properties - mass, size, mor- 
phology, gas content, star-formation rate, nuclear activity, 
characteristic velocity and metallicity. In particular, clus- 
tering studies connect galaxies to their unseen dark matter 
haloes and indicate how these were assembled. According 
to the current standard ACDM paradigm, galaxies form as 
gas cools and condenses at the centres of a hierarchically 
aggregating population of dark matter haloes, as originally 



outlined bv lWhite fc Reed (|l978l ). As smaller haloes fall into 
more massive ones, their central galaxies become satellites 
of these new hosts, occasionally merging at some later time 
into the new central galaxies in their cores. Thus, each halo 
contains a dominant galaxy at the bottom of its potential 
well, and a set of satellites which were the central galax- 
ies of smaller progenitors. Observational studies of satellite 
populations provide a check on this picture, indicating how 
central galaxy properties relate to halo mass, and how these 
properties are modified when a halo falls into a bigger sys- 
tem. 
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The abundance of satellites, their spatial distribution 
and their intrinsic properties are thus intimately bound up 
with halo merger histories, which are themselves closely re- 
lated to the underlying cosmology. For example, the evolu- 
tion of merger rates is sensitive to the cosmic matter den- 
sity, while the mass distribution of merging objects depends 
on the linear power spectrum of initial density fluctuations 
IjLacev fc Coldl 19931 ). Thus, satellite properties can, in prin- 
ciple, be used to constrain cosmological parameters. On the 
other hand, the physical processes driving galaxy evolution 
have strong effects on satellites. For example, their colours 
are affected by stripping of the gas reservoirs which supply 
star formation, and both gravitational and hydrodynamical 
processes can modify their structure, changing their mor- 
phology and partially or even totally disrupting them. Mod- 
ern evolutionary models for the galaxy population attempt 
to include such processes and can be tested by comparison 
with the abundances, colours and spatial distributions of 
satellites. 

The "missing satellite problem" is a particularly strik- 
ing example of how satellite galaxy studies can constrain 
cosmology and galaxy evolution. The problem highlights an 
apparent mismatch between the large number of self-bound 
subhaloes found in ACDM simulations of the formation of 
haloes like those of the Milky Way and M 31, and the much 
smaller number of satellite galaxies observed around these 
two Local Group galaxies JK lypin o t al.l Il999l : iMoore et al.l 
Il999l : iKravtsov. Gnedin. fc K lvpin 2004 '). This discrepancy 
has traditionally been addressed by claiming that photoion- 
isation and supernova feedback suppress cooling and star 
formation in low-mass haloes, so that only the most massive 
Milky Way subhal oes were able to make stars, the rest re- 
maining dark (e.g. Kauffmann. White, fc Guiderdonil 1992; 



Bullock. Kravtsov. fc Weinberg) l2000l: iBenson et allboO^ 



Somervilld l2002l : iMaccio et al.||2010l : Iguo et al.ll201lh . Re 
cent analyses of the kinematics of Galactic satellites sug- 
gest, however, that their dark matter haloes are not dense 
enough to corresp ond to the most massive subhaloes in 
a ACDM universe (IBoylan-Kolchin. Bullock, fc Kaplinghati 



|2012| : iFerrero et al.ll201lh . Such problems have led a num 
ber of authors to invoke warm dark matter (WDM) to 
eliminate subhaloes less massive than a few billion solar 
masses, and to re duce the central den s ity of subhaloes above 
this c ut-off (e.g. iMoore et all [JOOol: ISpergel fc Steinhar"dtl 



20001: lYos hida ct al.' '2000': 'Bode. O striker. fc Turokl l200ll : 
Zavala et al. 2009; LovcU ct al. 201^). Others claim that 
part of the problem may be in completeness of the 



observed satellite p o pulation (e.g., IWillman et al.l 2004: 



Simon fc Gehal l2007l: IKoposov et all l2008l : iToUerud et al] 



20081 : IWalsh. Willman. fc Jerienll200g ) 

Two kinds of methods are now commonly used to 
compare galaxy clustering in large redshift surveys to 
the predictions of high-resolution simulations of cosmic 
structure for mation. The Halo Occupa ti on Distribution 
(HOD ; e.g.. Ijing. Mo. fc Boernd Il998l: IJina^fc Boerner 



19981 : iPeacock fc SmithI I2OO0I: iMa fc Frvl I2OO0I: ISeliak 



200d: iBerlind and Weinberd I2OO2I : ICoorav fc ShethI |2002| : 
Zheng et al.l 120051 ) and t he closely related Con ditional Lu- 
minosity Function (CLF: lYang. Mo. fc van den Bosch 200j) 
approaches determine the central and satellite galaxy pop- 
ulations of haloes as a function of mass by optimizing the 
fit to abundance and clustering observations, typically lumi- 



nosity and correlation functions. In a non-parametric vari- 
ant, the observed abundance of central galaxies is matched 
directly to the simulated abundance of haloes to obtain 
a monotonic relation between galaxy luminosi ty an d halo 
mass (Abundance Matching, AM: iTas itsiom i et al.l 
Vale fc Ostrikej|2004l: IConrov. Wechsler. fc Kravtsovl 



2004 



2006 



Moster et al.ll2010l : iGuo et al.ll201Gl ). This relation can be 
used to populate both main and satellite subhaloes, if it is 
assumed to hold when satellites first fall into more massive 
systems. By construction, these methods fit observed lumi- 
nosity functions perfectly. The same is true for observed cor- 
relation functions for HOD and CLF, whereas these serve as 
a test for AM. None of these schemes explains why haloes of 
given mass have central galaxies with specific properties. 

In contrast, semi-analytic models use simplified rep- 
resentations of the relevant astrophysics to follow galaxy 
growth within the evolving dark halo population, and 
so attempt to predict the detailed pro perties of both 
centr a l and satellite galaxies (SAM; e.g ., White fc Frenk 
19911 : iKauffmann. White, fc Guiderdonil 1 19931: ICole et al 



1994: iKauffmann et al.1 Il999l: ISomerville fc Primackl I199I 



Springel et al.ll200ll : lKang et al.ll2005h . Here, adjustable pa- 
rameters correspond to the efficiencies of poorly under- 
stood processes like star formation or AGN feedback, so 
that the values derived from fitting to observation are in- 
teresting in their own right. In recent years, ever more 
detailed astrophysical models have been incorporated into 
ever larger and higher resolution simulations of dark matter 
evolution, leading to increasingly fait hful representation of 
the observe d galax y population (e.g. ISpringel et al.l l2005l: 



Croton et al. 2006: Bower et al] I2OO6I : iDe Lucia fc BlaizotI 



20071 : iGuo et al.ll20lil T 



Many studies of satellite galaxy populations have fo- 
cused on the Local Group (LG) because of the greater 
depth and detail with which n earby galax ies can be 
studied (see ICrebeJ (|2000l . l200ll . |2007| . I2OIII ) and refer- 
ences therein). Recent work has been particularly con- 
cerned with comparing ACDM predictions to the abun- 
dance and internal structure of dwarf spheroidal galaxies 
and to the orbital and internal prop erties of the Mag- 
ellanic Clouds (e.g. IKoposov et al.l 12009: Liu et al.l 2011 



Font et al. 201 
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e. g. IKoposov et al.l 1200 9: Li u et al.l 
ToUerud et al.l l201lUSales et al. 



|[ 



2011 



Bovlan-Kolchin. Bullock, fc Kaplinghaj [2oilh . Such stud- 
ies are limited by the fact that the LG contains only 
two large spirals, since considerable scatter is expected 
among the satellite popula t ions of similar m ass haloes (e.g. 
iBovlan-Kolchin et al.|[2010l : ICuo et al.ll201ll ). 

Beyond the LG, many observational studies of satel- 
lite populations have estimated their luminosity functions 
and their radial distribution around their primaries, but 
rather few have compared directly with theoretical expec- 
tations. Some studies have used redshift surveys to investi- 
gate the projected num b er density profiles of satellites (e.g. 
Vader fc Sandagel Il99ll : ISales fc Lamba^ l2005l : IChen et al.l 



20061 ) usually fitting power laws E(r) ~ r", and obtaining 



slopes a between —0.5 and —1.2. The availability of redshifts 
for all galaxies facilitates discrimination between satellites 
and background objects, but, for most objects, satellites are 
detectable only one or two magnitudes fainter than their 
primaries. An important application made possible by the 
redshifts is the measurement of mean dynamical masses for 
haloes as a function of central galaxy luminosity and colour 
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tearitsky et~aDll99 l 'l997^; IPrada et al1l2003l : IConrov et aD 
r200a . 120071 : iMore et al...201l1 '). The last of these papers finds 
that, at given luminosity, red central galaxies have more 
massive haloes than blue ones, but that this difference goes 
away if red and blue primaries are compared at the same 
stellar mass. We will return to this issue below. 

The abundance of satellites at magnitudes much fainter 
than their primaries is most easily studied by combin- 
ing a redshift survey with a photometric survey which 
catalogues galaxies several magnitudes below the spectro- 
scopic limit. In the absence of redshifts, it is not pos- 
sible, of course, to distinguish true satellites from back- 
ground galaxies. Results can therefore be obtained only 
by stacking large samples of primaries so that a statisti- 
cal substraction of the background population is possible 
llHolmberelll969l: iPhillipps fc Shanksfll987l: iLorrimer et al.l 



ll994l : ISmitli, Martinez, fc Graharnll2004l : lTal et al.ll2012h 71n 
particular, iLorrimer et al.l ( 1994 ) counted the number of 
faint images on Schmidt survey plates around primaries of 
known redshift, using a "bootstrap" method to remove the 
background and fitting the projected surface density to a 
power law ^{rp) ^ r^" , finding a ~ 0.9. They showed that 
satellites are more abundant and are concentrated to smaller 
radii around early-type primaries than around late-types. 

Most recent work has taken advantage of the enor- 
mous increase i n data provided b y the Sloan Digi t al Sky 



mous increase m gata proviaea D y tne aloan uigi t al oKy 
Survey (SDSS; I York et al.l |2000|). Iw einmann et al.l (|2006i ') 
used the group catalogue of lYang et al., (,i2007i ). constructed 
from the SDSS spectroscopic data, to study in detail how 
the properties of satellite galaxies depend on the colour, 
luminosity, and morphology of the central galaxy and on 
their inferred dark halo mass. They compared their obser- 
vational results_with a mock redshift survey based on the 
SAM of ICroton et al.l (|2006l ). finding significant discrepan- 
cies. In particular, the model overpredicted the number of 
faint satellites in massive haloes and produced too many 
red satellites. The fraction of b lue central galaxi e s was also 
too high at high luminosities. IWeinmann et al.l (|2006l ) ar- 
gued that the satellite problems most likely refiect an im- 
proper treatment of tidal stripping or of the truncation of 
star formation, while the central problem may reflect an 
overly simple treatmen t of dust or of AGN feedback. In 
IWeinmann et al.l (12011) a mock catalog based on the more 
recent model of lGuo et al.l (|201ll ) was compared with several 
nearby ga laxy clusters as well as with the group catalogue of 
lYang et al. , (,2 007j ) . Discrepancies were weaker than for the 
earlier model, but the predicted fraction of red dwarf satel- 
lites remains higher than in the Virgo cluster or in the group 
catalogue of lYang et al.l (|2007T ). although agreeing with the 
fractions found in the Coma and Perseus clusters. 

Studies of satellite galaxies based on both spec- 
troscopic and photometric data from the SDSS have 
been publis hed recent ly bv [Lares. Lambas. fc Dominguea 
ll201ll). bv IGuo et all (I2OIII) a nd bv iTal et all (|2012| )T 
iTal. Wake fc van DokkumI (|2012l ). The first of these investi- 



gated how the luminosity function and number density pro- 
file of satellites depend on the colour and luminosity of their 
central galaxy, finding the abundance of satellites to depend 
strongly on primary luminosity, and the faint-end slope of 
their luminosity function to be consistent wit h that of the 
field. Using similar datasets, IGuo et al.l (|201ll ) investigated 
the satellite luminosity function and its dependence on pri- 



mary luminosity, colour and concentration. Their satellite 
luminosity function estimates are not well fit by Schechter 
functions, tending to be flat at bright luminosities but very 
steep at fain t luminosities, apparently at odds with the 
conclusions of lLares. Lambas. fc Dominguea (|201H ). For the 
primary magnitude ra nge (My = —21. 25 ± 0.5) the mean 
luminosity function of IGuo et al.l (|201ll )is similar in shape 
to that of the MW and M31, bu t the abundance of satel- 
lites is about a factor two lower. iTal et al.l (|2012h studied 
satellites of SDSS Luminous Red Galaxies using, in partic- 
ular, the deep Stripe 82 data finding a luminosity function 
with a shallo w faint end slope and a very different shape 
from t hose of IGuo et all (|201ll ). lTaI. Wake, fc van DokkumI 
([2OI4) constructed radial number density proflles for these 
same s ystems, concluding that they are well fi tted by a NFW 
model (|Navarro. Frenk. fc White!! 19961 . 1 19971 ) on large scales 
while at small radii there is an excess of satellites compared 
with the NFW profile, which can be well dscribed by a Sersic 
model. Using data from th e Gala xy a n d Mass Assemb l y Sur - 
vey (GAMA: lD"river et al.ii200a . !20lII ). !Prescott et all|201l! ) 
also studied satellite number density profiles and red frac- 
tions as functions of projected separation and the masses of 
both satellite and primary, arguing that their results favour 
removal of gas reservoirs as the main m echanism quench- 
ing st ar formation in satellites. Finally, iNierenberg et al.! 
(|2012h use HST data from the Cosmological Evolution Sur- 
vey (COSMOS) to study similar issues for smaller samples 
of galaxies, but out to z ~ 0.8. 



In the present paper we return to many of these ques- 
tions, using the full SDSS spectroscopic and photometric 
databases i n conjunct i on wit h the galaxy population sim- 
ulations of !Guo et al.! (|201l! , hereafter Gil). The simula- 
tions allow us to compare expectations based on our cur- 
rent understanding of galaxy formation in a ACDM universe 
with the observed dependences of satellite luminosity, stel- 
lar mass and colour on primary galaxy properties. By using 
mock catalogues from the simulations, we are able to explore 
how satellite populations relate to the dark matter haloes in 
which they are embedded, to gain insight into the effect of 
physical processes like quenching and tidal stripping on their 
observable properties, and to explore how the observational 
criteria defining isolated primary galaxies impact the clus - 
tering of other galaxies around them (see iFall et al.! (|l97d ) 
for an old example of the potential strength of such effects). 



We describe the datasets we use and the selection cri- 
teria which define our primary and satellite galaxy samples 
in section [2] In section |3l we introduce our background sub- 
traction method. We present our SDSS results and compare 
them directly with the Gil simulation in sections |4] and [5| 
Further discussion and comparison with previous work is 
given in our concluding section. An appendix describes a 
variety of tests for systematics in the SDSS photometric 
data and in the techniques we use to correct satellite counts 
for contamination by foreground and background galaxies. 
Throughout this paper, we convert observational to intrinsic 
properties assuming a cosmology with Qm = 0.25, Q,a ~ 0.75 
and h — 0.73. We quote all masses in units of A/q rather 
than h~'^MQ. 
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2 DATA AND SAMPLE SELECTION 

2.1 Primary Selection 

We wish to study the satellite populations of bright iso- 
lated galaxies out to distances ~ 0.5 Mpc. We begin by 
considering all galaxies brighter than r = 16.6 (r-band ex- 
tinction corrected Petrosian magnitude) in the spectroscopic 
galaxy catalogue of the New Yor k University Value Added 
Galaxy Catalog fNY U-VAGCfl (lAihara et al.ll2011^ . This 
catalogue was built bv lBlanton et al.l (|2005l ) on the basis of 
the seventh Data Release of the Sl oan Digital Sky Survey 
('SDSS/DR7: lAbazaiian et alllJOOgh . This apparent magni- 
tude limit provides us with a parent catalogue of 145070 
objects. We select isolated galaxies from this sample by re- 
quiring that there should be no companion in the spectro- 
scopic sample at r^ < 0.5 Mpc and |A2| < 1000 km/s that 
is less than a magnitude fainter in r than the central object, 
and no companion at rp < 1 Mpc and \Az\ < 1000 km/s 
that is brighter than it. These isolation criteria reduce our 
sample to 66285 objects. 

The SDSS spectroscopic sample is incomplete, because 
observing efficiency constraints made it impossible to put a 
fibre on all candidates or to re-observe objects where an 
initial spectrum yielded an unreliable redshift. The com- 
pleteness varies with position on the sky and has a mean 
of 91.5% for our parent sample. Thus ~ 10% of our "iso- 
lated" galaxies will not, in fact, be isolated according to our 
criteria, because their companion was missed by the spec- 
troscopic survey. To eliminate such systems we apply an 
additional cut using the SDSS photome tric data. The pho - 
tometric redshift 2 catalogue (photoz2; ICunha et al.l [2009|) 
on the SDSS website provides redshift probability distribu- 
tions for all galaxies in the SDSS footprint down to apparent 
magnitude limits much fainter than we require. These dis- 
tributions are tabulated for 100 redshift bins, centered from 
zi = 0.03 to Z2 = 1.47 with spacing dz = 1.44/99. We find 
all the objects in our candidate isolated galaxy list which 
have a companion in the photoz2 catalogue satisfying the 
above projected separation and magnitude difference crite- 
ria, and we discard those where the companion has a pho- 
tometrically estimated redshift distribution compatible with 
the spectroscopic redshift of the primary. Our definition of 
"compatible" is that the probability for the companion to 
have a redshift equal to or less than that of the primary ex- 
ceeds 0.1. Apparent companions which fail this test usually 
do so because their colours are too red to be consistent with 
a redshift as low as that of the primary. After applying this 
additional cut, 41883 objects remain in our isolated galaxy 
sample. 

Finally, we exclude any object for which more than 20% 
of a surrounding disc of radius Vp — IMpc lies outside the 
survey footprint. Such objects could have bright companions 
which are not included in the SDSS databases. To evaluate 
these overlaps we made use of the set of "spherical poly- 
gons" provided on the NYU-VAGC website. These account 
both for the survey boundary and for masked areas around 
bright stars. We generate a large number of points uniformly 
and randomly over the 1 Mpc disc surrounding each galaxy 
and discard any which lie outside the .survey boundary. A 



http://sdss.pliysics.nyu.edu/vagc/ 



galaxy is eliminated from the sample if more than 20% of its 
points are discarded in this way. This last cut removes about 
1.5% of our objects, leaving a final sample of 41271 bright 
isolated galaxies. The set of randomly generated points sur- 
rounding each of these "primaries" is kept for later use when 
estimating background corrections to the counts of its faint 
companions (see below). 

Figure [l] compares the distributions of colour, Petrosian 
half-light radius -Rso, concentration C — Rgo/Rso and stellar 
surface mass density fi-,, = Mt/2-KB3,Q for our parent galaxy 
sample (all 145070 galaxies with r < 16.6) and for our 41271 
isolated primaries, separated into bins of galaxy stellar mass. 
These quantities were taken directly from the NYU-VAGC 
catalogue. The stellar masses were estimated by fitting stel- 
lar population syn thesis models to the K-corrected galaxy 
co lours assuming alChabrien (|2003h initial mass function as 
in iBlanton fc Roweij ( 20071 ) . The sensitivity of the stellar 
masses to assumptions underl ying the estimatio n technique 
is explored in the Appendix of lLi fc White! ([20091). Each stel- 
lar mass bin is a factor of two wide, and we show data for 
the five mass bins on which we will concentrate our analy- 
sis throughout the rest of this paper. The numbers in the 
lower right of the -R50 plots indicate the number of isolated 
primaries in each mass bin. Volume corrections have been 
applied to all the histograms in this figure by calculating 
the total volume Vmax of the survey over which each indi- 
vidual galaxy would be brighter than the flux limit, r = 16.6, 
and accumulating counts weighted by 1/1/max. It turns out 
that our isolated primaries are slightly bluer than the par- 
ent sample, particularly at low masses. In addition, they are 
slightly more concentrated than the parent sample. Our se- 
lection procedure appears to have no significant effect on 
the distributions of the other two properties, and the same 
is true for the redshift distributions which we do not show. 
To a very good approximation our isolated primary galax- 
ies are typical objects of their stellar mass, although as we 
will see below, our selection has a strong infiuence on their 
relation to their environment. Vertical dashed lines in the 
colour plots indicate the split we adopt when separating our 
primaries into red and blue populations. This split is slightly 
dependent on stellar mass. 

2.2 The mock catalogue 

In the analysis and interpretation of our results for satellite 
galaxy populations we will make considerable use of mock 
catalogues built from the galaxy formation simulations of 
IGuo et al.l (|201ll . hereafter Gil). These are implemented 
on two very larg e dark matter simula tions, the Millennium 
Simulation (MS; ISpringel et al.1 120051') and the M illennium- 
II Simulation fMS-II: iBovlan-Kolchin et al.ll2009h . The MS 
follows the evolution of structure within a cube of side 
500/i~^Mpc (comoving) and its merger trees are complete for 
subhaloes above a mass resolution limit of 1.7 x 10^''/i~^Mq. 
The MS-II follows a cube of side 100/i~^Mpc but with 
125 times better mass resolution (subhalo masses greater 
than 1.4 X 10**/i"^Mq). Both adopt the sa me WMAPl- 
based ACDM cosmology (jSpergel et al.l 120031 ) with param- 
eters h = 0.73, flm = 0.25, D.A = 0.75, n = 1 and o-g = 0.9. 
These are outside the region preferred by more recent anal- 
yses (in particular, ag appears too high) but this is of no 
consequence for the issues we study in this paper. For consis- 
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Figure 1. Volume weighted distributions for the parent sample of SDSS/DR7 galaxies with r < 16.6 (black curves) and for our sample 
of isolated galaxies (red curves). Different panels within each set refer to different ranges of log M*/M0 as labelled. Each set corresponds 
to a different property. Top left: "'^(g — r) colour - the vertical dashed line in each panel shows the colour separating red and blue 
populations. Top right: Concentration C = Rgo/fiso- Bottom left: Petrosian half-light radius R50 (iii kpc) — the numbers at bottom 
right of each panel indicate the numbers of galaxies in the parent (black) and isolated samples (red) shown. Bottom right: Stellar 
surface mass density /i* = Mt/27rij|Q. 



tency, we will adopt this cosmology when quoting numbers 
in the rest of this paper. 

In Gil's galaxy formation model, the uncertain star for- 
mation and feedback efficiencies are tuned to produce close 
fits to the stellar mass, luminosity and autocorrelation func- 
tions of low redshift galaxies as inferred from the SDSS. To- 
gether with their high resolution (particularly for the MS-II) 
and large volume (for the MS) this makes them ideal for our 
purposes in this paper. Here we use the publicly available 
data from http://www.mpa-garching.mpg.dG/millennium 
We project the simulation boxes in three orthogonal direc- 
tions parallel to their x, y and z axes. In each projection 



we can assign each galaxy a redshift based on its "line-of- 
sight" distance and peculiar velocity, and we can select iso- 
lated primaries using criteria which are directly analogous 
to those used for the SDSS (though we do not need to worry 
about the complications due to completeness and boundary 
issues). In addition to observables like luminosities, colours, 
sizes and morphologies, the simulation databases provide in- 
formation which is not directly accessible for real galaxies 
(e.g. halo mass, environment type and full 3D position and 
peculiar velocity). These can give insight into the nature of 
the isolated primary sample we have selected. 

All the SDSS luminosities and colours we use in this pa- 
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per are rest-frame quantities K-corrected to the ' r band. 
The absolute magnitudes in our mock catalogues are in the 
true z = r band, because this is what is directly pro- 
vided by the database and the difTerence between the '^r 
and r bands is too small to significantly afi'ect luminosities. 
We do, however, transform the database {g — r) colours to 
the °'^(.9 — r) band usin g the empirical fitting formula of 
iBlanton fc Roweij (|2007l ) because here the shifts seem large 
enough to cause (minor) differences. 

Structure in the MS and MS-II is characterized us- 
ing Friends-of Friends (FoF) groups partitioned into sets 
of disjoint self-bound subhaloes. The subhalo populations 
at neighboring output times are linked to build merger trees 
which record the assembly history of all nonlinear structures 
and provide the framework for the galaxy formation simu- 
lations. In these simulations galaxy evolution is affected by 
environment in several ways. The galaxy at the centre of 
the most massive subhalo of each FoF group (which usually 
contains most of its mass) is considered the "central galaxy" 
and is the only one to accrete material from the difi^use gas 
associated with the group. When evolution joins two FoF 
groups, Gil continue to treat the galaxy at the centre of the 
less massive subhalo as a central galaxy until it falls within 
the nominal virial radiu^J of the new FoF group. After this 
point the infalling galaxy is considered as a "satellite", the 
mass of its dark halo starts dropping as a result of tidal 
stripping, and its diff'use gas is assumed to be removed in 
proportion to the subhalo dark matter. Such satellites may 
later lose their subhaloes entirely through tidal disruption. 
At this point they are either disrupted themselves or (more 
commonly) they become "orphan satellites" which continue 
to orbit until dynamical friction causes a merger with their 
central galaxy. In this paper we will follow Gil, considering 
together the two kinds of satellites (with and without a dark 
matter subhalo) and the two kinds of centrals (in dominant 
and in newly accreted, distant subhaloes). 

The left panels of figure [2] compare halo mass distri- 
butions as a function of stellar mass for isolated simula- 
tion galaxies to those of their parent population. Here, halo 
mass is M200 of the FoF group for the dominant central 
galaxy and all the satellites, and Minf for the other "cen- 
trals", where Mi^f is the Af200 of the old FoF group of a 
newly accreted central just prior to infall. Our isolation cri- 
teria have a substantial effect on these mass distributions, 
eliminating a high-mass tail which is particularly evident 
for lower stellar mass galaxies. This tail is due to the satel- 
lites which, as we will show more explicitly below, are very 
effectively excluded by our criteria. 

The right panels of figure[2]compare colour distributions 
for these same simulated galaxy populations. Just as for the 
real SDSS galaxies (figure [ij our isolation criteria bias the 
distributions to bluer galaxies because central galaxies have 
ongoing gas accretion and so are more actively star-forming 
than satellites. Selection induces larger shifts for the sim- 
ulated distributions than for the real ones, however, and 
the bimodal nature of the colour distributions is more ob- 



^ We define this as r2oo the radius of the sphere centred on the 
gravitational potential minimum of the FoF group within which 
the mean density is 200 times the critical value. 4^200 is then the 
mass within this sphere. 



vious in the sim u lation . As already discussed by Gil and 
IWeinmann et al.l (|201ll ) this refiects the facts that the red 
and blue sequences are more sharply defined in the simula- 
tion than in reality and that satellite galaxies appear to be 
too uniformly red. Note also that the red and blue popula- 
tions separate at bluer colours (indicated by dashed vertical 
lines) in the simulation than in the SDSS data, particularly 
at high mass. This appears to be a consequence of the stellar 
population synthesis models used in the simulation, together 
with the fact that stellar metallicities are too low for high- 
mass galaxies (see iHenrigues fc Thomasll2010l ). 

The left panels of figure O illustrate how well our iso- 
lation criteria select central galaxies, at least in the sim- 
ulation. The fraction of centrals according to the defini- 
tions of Gil is plotted as a function of stellar mass in 
the lower panel, with the black curve referring to all iso- 
lated galaxies and the red and blue curves referring to 
the red and blue subpopulations. For the isolated popu- 
lation as a whole and for its red subpopulation, the con- 
tamination by satellites maximizes at just over 2%. Slightly 
larger contamination occurs at high mass in the blue sub- 
population, but such massive blue galaxies are in any case 
very rare (see below). In contrast, the curves in the upper 
panel show that only 65% of the parent population in the 
Gil model are centrals at log M^/M© = 10, about 80% 
are centrals at logAf,t/M0 = 11 and 88% are centrals at 
logM*/Af0 = 11.5. For blue galaxies the central fraction is 
above 90% at all masses, while for logAf*/Af0 < 10.3 most 
red galaxies are satellites. The application of our isolation 
criteria to the simulated galaxy population of Gil thus re- 
sults in an extremely pure sample of central galaxies. 

We cannot be sure, of course, that the elimination of 
satellite galaxies is as effective in the SDSS samples as in 
the simulation. We can, however, check that the separation 
into red and blue subpopulations matches as a function of 
stellar mass, both for isolated galaxies and for their parent 
population. This comparison is shown in the right panel of 
figure |3] Here red curves indicate the red fraction as a func- 
tion of stellar mass for isolated galaxies, while black curves 
indicate the same quantity for the parent sample. In each 
case, solid curves are the observational result from SDSS 
and dashed curves are for the simulation of Gil. For the 
observations we have again used 1/Knax weighting to en- 
sure that the plotted quantity is appropriate for a volume- 
limited sample. For the primary stellar masses relevant for 
this paper ( log Af^/Af© > 10.2) the agreement between ob- 
servation and simulation is almost perfect. At lower stellar 
masses there are too many red galaxies in the simulation, 
again refl ecting the problem not ed above and discussed in 
detail bv IWeinmann et al.l l|201ll l: at low masses the simu- 
lated satellite galaxies are too uniformly red. 

Another comparison of the effects of our isolation crite- 
ria on the observed and simulated galaxy samples is shown in 
Figure 2] The thick black solid line here indicates, as a func- 
tion of stellar mass, the fraction of galaxies from the parent 
sample which remain in our final sample of isolated SDSS 
galaxies. The thick black dashed line shows the result when 
analogous isolation criteria are applied to our Gil mock 
catalogues. The agreement of observation and simulation is 
again quite good, though not as perfect as in the right panel 
of figure [3] Note, however, that perfect agreement is not 



© 2011 RAS, MNRAS 000, [T]-?? 



o 0.1 - 



o 0.1 



2 




111. 4-11. 7 



LIO. 8-11.1 



Satellite abundances 7 



JLll.l-11.4 



JLlO.5-10. 



.10.2-10.5 



_all 
jsolated 



11 12 13 14 15 11 12 13 14 15 

loglo(Mhalo/Mo) loglo(Mhalo/Mo) 



0.5 
g-r 



0.5 



Figure 2. Left: Halo mass distributions for isolated galaxies (red curves) in five disjoint stellar mass bins and for the parent populations 
in the simulations of Gil from which they were drawn (black curves). Right: Colour distributions for these same sets of simulated 
galaxies. Vertical dashed lines indicate the colour at which we separate red and blue populations. 
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Figure 3. Left: The fraction of isolated galaxies (bottom panel) and of their parent sample of all galaxies (top panel) in the galaxy 
formation simulation of Gil which are classified as centrals rather than satellites according to the Gil criteria which we also adopt here 
(see the text for details). Black, red and blue lines indicate the fractions as a function of stellar mass for all, for red and for blue galaxies 
respectively. Right : The fractions of galaxies which are red (rather than blue) are shown as a function of stellar mass for isolated galaxies 
(red curves) and for the parent sample of all galaxies (black curves). The solid curves in each case represent observed results for the 
SDSS (after weighting to represent values for volume-limited samples) while the dashed curves are for the galaxies in the simulation of 
Gil, Note the excellent agreement for logM*/MQ > 10.2 the range of interest for this paper. 



expected, since the underlying SDSS sample is magnitude- 
limited whereas the Gil sample is volume- limited. 

This same plot provides a convenient way to summaries 
the relative effects of the spectroscopic and photometric cat- 
alogues in defining our observed sample of isolated galax- 
ies. The thin red solid line shows the fraction of objects 



which are isolated relative to the spectroscopic sample, but 
not necessarily relative to the photometric sample. At most 
masses, use of the photometric catalogue increases the num- 
ber of objects with an identified "companion" by 30 to 50%. 
This has a relatively small effect on the isolated fraction at 
high mass, since the great majority of massive galaxies are 
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isolated by our criteria. At log M^/M© = 10.3 it reduces 
the number of isolated galaxies by almost a factor of two, 
however, because the majority of such SDSS galaxies have 
a companion by these same criteria. These numbers suggest 
(and we have checked in more detail) that our photometric 
rejection step is conservative, in that a significant fraction 
of the photometric "companions" (about one third) are, in 
fact, galax;ies at a significantly different redshift which are 
projected on top of the primary. This does not matter for 
the analysis of this paper - it just causes us to end up with a 
slightly smaller sample of isolated primaries than if we had 
been less conservative. 

Finally we can use the simulated catalogue to see how 
our isolation criteria affect the numbers of satellite and cen- 
tral galaxies in our samples. The green and blue dashed 
lines in Figure |4] show as a function of stellar mass the frac- 
tions of centrals and satellites which pass our isolation crite- 
ria. At high mass, most centrals are indeed isolated, but for 
log Af*/Af0 < 10.6 more than half of them are rejected be- 
cause of apparent companions. In most cases these compan- 
ions are actually the central galaxies of other haloes. On the 
other hand, for log M^/M© < 11 almost all satellite galaxies 
are found to have companions according to our criteria. With 
increasing mass the number of apparently isolated satellites 
grows, reaching 30% at the highest mass. This is because 
these objects lie in very massive galaxy clusters and so can 
be projected more than 1 Mpc from their associated cen- 
tral galaxies. Note, however, that the actual number of such 
nicissive satellites is very small (see figure |3]). It is the inter- 
play of these different effects which produces the very high 
purity (i.e. central galaxy fraction) at all masses in our final 
isolated galaxy sample. 




11 11.5 

logi„(M./Mj 

Figure 4. The fraction of galaxies that are selected by our iso- 
lation criteria as a function of stellar mass. The thick black solid 
line gives the result for our final sample of isolated SDSS galax- 
ies, while the thick black dashed line is the corresponding result 
for simulated galaxies in the Gil model. These agree moderately 
well over the full mass range shown. The thin solid red line shows 
the fraction of retained SDSS galaxies after demanding isolation 
relative to the spectroscopic sample but before additionally re- 
quiring isolation relative to the photometric sample. Green and 
blue dashed curves indicate the fractions of central and satellite 
galaxies which remain after applying our isolation criteria to the 
Gil simulated galaxy sample. 



2.3 The photometric catalogue for satellites 

Both the spectroscopic catalogue and the photoz2 photo- 
metric catalogue used in the last section are based on DR7, 
the seventh release of SDSS data. When identifying and ver- 
ifying the isolation of our isolated primary galaxies, only ob- 
jects with apparent magnitude brighter than r — 17.6 had 
to be considered. This is far above the magnitude limit of 
the SDSS photometry, and the accuracy of the DR7 mag- 
nitudes/colours and of the derived photo-z distributions is 
well tested at these magnitudes. When compiling counts 
of faint satellite galaxies around these primaries, we need 
to go to the SDSS limit for reliable photometry, however, 
which we take to be r = 21. At this limit, improvements 
to the photometry pipeline have continually enhanced its 
reliability, identifying (and correcting when possible) sys- 
tematic artifacts which can have a significant infiuence on 
our analysis. We therefore use the photometric catalogue 
from SDSS/DR^ (JAihara et al.ll201ll ) when compiling satel- 
lite counts. To be specific we created a reference photo- 
metric catalogue by downloading objects that are classi- 
fied as galaxies in the survey's primary object list, and 
that do not have any of the flags BRIGHT, SATURATED, 
SATUR.CENTER or NOPETRO_BIG set. This follows the 
selection criteria for the DR7 photoz2 catalogue used above, 
It is important to note that the magnitude we use for 



^ http://skyscrvice.pha.jhu.edu/casjobs/ 



satellites is the so-called SDSS model magnitude, which, at 
faint apparent magnitudes, is claimed to have the highest 
signal-to- noise of all the alternatives catalogued. In contrast, 
when defining our primary sample and checking its isola- 
tion we used Petrosian magnitudes as listed on the NYU- 
VAGC website. We continue to use these magnitudes for 
the primaries below. When we quote colours, these are al- 
ways measured within a well-defined aperture related to the 
Petrosian radius, and are rest-frame quantities K-corrected 
to the "'^{g — r) band for both primaries and satellites. All 
absolute magnitudes quoted for SDSS galaxies are also in 
this same rest-frame band. 

We have carried out a variety of tests for systematics in 
the SDSS photometry, checking completeness and the qual- 
ity of star-galaxy separation by comparing with much deeper 
HST data, and quantifying systematic biases in the magni- 
tudes of faint images (e.g. the satellites) in the neighbor- 
hood of substantially brighter images (e.g. the primaries). 
Detailed descriptions of these tests are given in an Appendix, 
along with details of tests of our procedures (outlined in the 
next section) for counting apparent companions around our 
primary galaxies and correcting for the (usually dominant) 
contribution from unrelated foreground and background ob- 
jects. 
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3 SATELLITE COUNTING METHODOLOGY 

We want to study the abundance of satellites as a func- 
tion of their luminosity, stellar mass and colour around our 
sample of isolated bright primaries, and to see how these 
abundances depend on the stellar mass and colour of the 
primary. We will need to use the SDSS data down to their 
reliable photometric limit (which we take to be r = 21) and 
as a result the great majority of the potential satellites do 
not have a spectroscopically measured redshift. Hence it is 
necessary to count all apparent neighbours around our pri- 
maries and to correct statistically for unassociated objects 
which happen to be projected near them. The number of 
such (mainly) background objects can substantially exceed 
the number of true satellites, so it is important to take con- 
siderable care in making these corrections. We adopt the 
following procedure. 

For each isolated bright galaxy, we identify all photo- 
metric galaxies with apparent projected separation rp < 
0.5 Mpc and we accumulate counts in bins of projected sep- 
aration rp, apparent magnitude r and observed colour g — r. 
From the count in each {rp,r,g — r) bin, we subtract the 
expected number of background galaxies which we take to 
be N{r,g — r)A{rp,z)f/Atot, where N{r,g — r) is the total 
number of galaxies in the {r,g — r) bin in the full photomet- 
ric catalogue, Atot is the solid angle of the survey footprint, 
A{rp, Zpri) is the solid angle corresponding to the annular rp 
bin at the redshift Zpri of the primary galaxy, and / is the 
incompleteness factor, the fraction of this annulus which lies 
within the survey footprint (we estimate this using the ran- 
dom points generated around the position of each primary 
during the selection process - see above). We then use the 
redshift of the primary galaxy to convert observed appar- 
ent magnitudes and colours into rest-frame luminosities and 
colours (in the "'^r and ^'^{g — r) systemjj and we transfer 
the background-subtracted satellite counts from our narrow 
bins of r and g — r into substantially broader bins of the rest- 
frame quantities. Finally we average these counts for each 
rp bin over the set of all primaries in the desired range of 
stellar mass (and sometimes colour) and we sum the result 
over the desired range in rp. Uncertainties in the resulting 
numbers are estimated from the scatter among results for 
100 bootstrap resamplings of the set of primaries. 

Some apparent companions are too red to be at the red- 
shift of the primary galaxy. It is useful to exclude them when 
accumulating counts since they add noise without adding 
signal. Hence, we exclude all bins redder than ^'^{g — r) = 
0.0321ogjQ A'f* -f 0.73, a fit to the upper envelope of the distri- 
bution of rest-frame colour against stellar mass for galaxies 
of measured redshift. The stellar mass Nh of the apparent 
companion is estimated by assuming it to be at the primary's 
redshift and adopting 

-^r)^ + 4.1183''-^(p-r)-0.7837 (1) 



{M/L)r = -1.0819" 



This empirical relation is a fit to a flux-limited 
(r < 17.6) galaxy sample from the NYU-VAGC web- 
site for which stellar masses were estimated from the K- 
corrected galaxy colours by fitting stellar population synthe- 



* We use the empirical fitting formula of IWestra et al.l 1120101 ) 
which gives the K-correction as a function of redshift and observed 
colour. 



sis models assuming a | Chabrieri (|2003h initial mass function 
(jBlanton fc Roweis^l2007^ . For this sample the 1-a scatter in 
{M/L)r of this simple relation is about 0.1. 

The photometric catalogue we use is complete down to 
an r-band apparent model magnitude of 21. This limit cor- 
responds, of course, to different satellite luminosities and 
stellar masses for different primary redshifts and different 
satellite colours. In order to ensure that our samples are 
complete when compiling satellite luminosity functions, we 
allow a particular primary to contribute counts to a partic- 
ular luminosity bin only if the K-corrected absolute lumi- 
nosity corresponding to r = 21 for a galaxy at the redshift 
of the primary and lying on the red envelope of the intrinsic 
colour distribution is fainter than the lower luminosity limit 
of the bin. Thus only the nearest primaries will contribute to 
the faintest luminosity bins of our satellite luminosity func- 
tions, and different numbers of primaries will contribute to 
each bin. We follow an exactly analogous procedure when 
compiling stellar mass functions for satellites. 

For each individual satellite luminosity/stellar mass bin, 
this treatment is equivalent to imposing an upper limit on 
primary redshift. Thus there is a maximum volume which is 
surveyed for satellites in the jth luminosity/stellar mass bin 
which we denote VmaxMn.j- 

On the other hand, our primary sample is fiux-limited 
at r = 16.6, so for brighter satellite bins where VmaxMnj 
is large, intrinsically faint primaries will not be visible to 
the redshift limit. The effective volume surveyed is then 
Vmax.pri.i, the total survey volume over which the ith pri- 
mary would lie above the flux limit. Note that because of 
K-corrections, this volume depends on the intrinsic colour 
of the primary as well as its intrinsic luminosity. We present 
our final results in the form of the mean number of satellites 
per primary for a volume-limited primary sample 



Nsat,i — 



^i-^^ sat,i,j / Vm 



sa/F™ 



where 



K. 



i[K. 



,., V, 



max, bin, j \ 



(2) 



(3) 



Thus, for satellite luminosity/stellar mass bin j, we sum 
satellite counts over all primaries i that are within VmaxMn.j- 
At the bright end of the satellite luminosity function, we 
expect Vmax,ij ~ Vmax,pri,i becausc Satellites are less than 
4.4 magnitudes fainter than their primaries and so can be 
seen around all primaries. For intrinsically faint satellites, 



however, Vm 



Vr 



Tnax,hin,j 



because primaries can be seen 



to well beyond the distance at which r = 21 for the satellites. 
We apply similar selection criteria to our mock cata- 
logues based on Gil. Since we know the absolute magni- 
tude, rest-frame colour and stellar mass of all galaxies, the 
background subtraction can be carried out directly using 
rest-frame quantities. In addition, the effective depth of the 
satellite catalogue is the same for all primaries so that all 
primaries can contribute to all satellite luminosity or stellar 
mass bins and we do not need any weighting in order to ob- 
tain proper volume-weighted statistics. Note that we do not 
use any information about the redshift difference between 
primary and apparent companion, so projection effects oc- 
cur over Az = 50, 000 km/s corresponding to the side of the 
Millennium Simulation "box". Both for the simulation and 
for the real SDSS data we use a global background estimate 
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based on the full survey in order to minimize the statistical 
uncertainty in the correction. 

Previous work often est imated a background density lo- 
cally for each primary (e.g. iLorrimer et al.lll994l : IGuo et al.l 
[2OIII). This not only substantially increases the noise, it also 
introduces a significant bias since galaxies are correlated on 
all scales and the "background" fields are, in fact, expected 
to have faint galaxy densities significantly above the mean. 
The extent of the bias is, in fact, strongly dependent on the 
isol ation criteria an d can be of either sign (see, for e-y^ajn- 
ple. lFall et aLlll976l ). The large-scale uniformity of the SDSS 
photometry is such that there appears to be no advantage to 
adopt a local estimate, provided all photometric quantities 
are properly corrected for Galactic extinction. The accuracy 
of our method is confirmed by the accurate convergence to 
unity at large angular scale in the left panel of figure lAll 
and by a variety of tests which we describe in detail in the 
Appendix. 



4 LUMINOSITY AND MASS FUNCTIONS OF 
SATELLITE GALAXIES 

4.1 Luminosity functions 

In figure [S] we present "'^r-band luminosity functions for 
satellites projected within 300 kpc of their primaries, except 
for the faintest bin, where the halo virial radius is much 
smaller than 300 kpc and we estimate the luminosity func- 
tion within 170 kpc in order to increase the signal-to-noise. 
In the left panel, data points connected by solid lines show 
our observational results for SDSS/DR8. As indicated in the 
legend, colours encode the range in logM+/Af0 of the pri- 
maries contributing to each luminosity function estimate. 
Error bars are derived by bootstrap resampling the primary 
sample. At the faint end we lose higher redshift primaries 
because of the apparent magnitude limit of our photomet- 
ric catalogue. We do not plot data for bins with fewer than 
eight primaries. It is evident that satellite numbers increase 

strongly with primary mass. 

The black solid line in the left panel shows a lSchechtej 

(|l976r ) function fit to the data for the most massive pri- 
maries. We have fixed the characteristic luminosity and 
faint-end slo pe to be those of the S DSS field luminosity func- 
tion in "V (|Montero-Dorta fc Pra da 2009). The resuh is a 
moderately good fit to the satellite data, although these ap- 
pear steeper than the field luminosity function at the faintest 
magnitudes. The satellite luminosity functions clearly be- 
come steeper for lower mass primaries. Power-law fits to 
the faint-end data are shown as dashed lines in the figure 
and give the slopes a listed in the first row of table [T] The 
bright end of each function is cut off by our requirement 
that every satellite be at least one magnitude fainter than 
its primary. We therefore exclude one or two of the brightest 
points when making these fits. Specifically, we find the me- 
dian absolute magnitude Mmed for primaries in each mass 
range, and we include only points for which the correspond- 
ing absolute magnitude bin lies entirely below Mmed-l-l. The 
ranges fitted in each case are indicated by the extents of the 
dashed black straight lines|j The faint-end slope decreases 



We have checked that the remaining points are indeed unaf- 



from a ~ -1.2 for logM*.p/Af0 ~ 11.5 to a ~ -1.6 for 
log Mv^.p/Mq ^ 10.3. 

IGuo et all (|201ll ) divided their primaries into three 
luminosity ranges {Mr = -23.0 ± 0.5, -23.0 ± 0.5 and 
—23.0 ± 0.5) and compiled satellite luminosity functions in 
bins of satellite-primary magnitude difference rather than 
satellite absolute magnitude. They quote faint-end slopes of 
-1.45, -1.725 and -1.96 for these three sets of primaries, with 
the fainter primaries having steeper luminosity functions. 
For the largest measured magnitude differences (Am ~ 8) 
they found similar numbers of satellites independent of pri- 
mary luminosity. These a values are substantially more neg- 
ative than ours, particularly for the faintest primaries. In or- 
der to compare with their results, we adopt similar isolation 
criteria, we take the same ranges of primary luminosity, and 
we also accumulate satellite number as a function of mag- 
nitude difference. Fitting the faint-end sl ope over the sam e 
satellite magnitude range as in figure 7 of lGuo et al.1 (|201H ). 
we find a values of -1.1 89, -1.376 and -1.5 88, substantially 
shallower than those of IGuo et al.l (|201ll ) and quite com- 
patible with those we quote in Table [T] Detailed tests show 
this inconsistency to be due part l y to t he local background 
subtraction scheme of lGuo et al.l ( 20 111 ) which removes part 
of the signal q but mainly to the fact that they use model 
magnitudes K-corrected to z = for their primaries, rather 
than the "'^r Petrosian magnitudes which we use here. 

In the right panel of figure [5] the small symbols 
joined by solid lines show analogous results for the galaxy 
formation model of Gil, based on the Millennium and 
Millennium-II simulations. Points brighter than Mr = — 18 
are MS data. At fainter magnitudes, resolution effects cause 
the MS to underestimate galaxy abundances and we take our 
data from the MS-II. (The two simulations agree very well in 
the range — 18 > Mr > —20.) To facilitate comparison, we 
replot the SDSS data from the left panel as filled triangles. 
Agreement of model and observation is fair but far from 
perfect. The simulation overpredicts the number of satel- 
lites around the most massive primaries by 25 to 50% for 
Mr > —19.5. In the two lower primary stellar mass ranges, 
the simulation underpredicts the number of satellites by 20 
to 30% . As we will see below, the latter discrepancy reflects 
a problem with the colours of the simulated satellites rather 
than with their stellar masses. The black dashed lines in this 
panel are the "held" luminosity function for the full simula- 
tions renormalized to fit the satellite data for each primary 
mass range. Here also there is a trend for the faint-end slope 
to be steeper for satellites than in the field, but the effect is 
much less marked than for the SDSS data. Furthermore the 



fected by rebinning our data as a function of the r-band absolute 
magnitude of the primaries. In this case we know exactly which 
bins are unaffected by our isolation criterion. When primary abso- 
lute magnitude and stellar mass are matched appropriately, the 
resulting sateUite luminosity functions match those of figure \5\ 
very closely over the full range used to determine the faint-end 
sl ope and are unaff ected by the isolation criterion over this range. 
6 IGuo et al.l 1I2OIII) used photometry from SDSS/DR7 while our 
own tests, similar to those discussed in section Al of the Ap- 
pendix, showed to suffer from substantially more serious system- 
atics for faint images close to brighter ones than is the case for 
the SDSS/DR8 catalogues used here. 
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Figure 5. Luminosity functions in the "■^■r band for satellites of primary galaxies in five disjoint ranges of logM*/MQ, as indicated 
in the legend. Satellites are counted within a projected radius of 300 kpc, except for the lowest mass range where we count within 
170 kpc. In the left panel, the data points connected by solid lines give observ ational results fo r SDSS/DR8 with error bars estimated 
by bootstrap resampling of the primary sample. The solid black lino is a fit of a lSchechterl 1119761 ) function to the most massive bin, with 
the characteristic luminosity and faint-end slope constrained to match those of the SDSS field luminosity function. Dashed black lines 
are power-law fits excluding one or two of the brightest points for each mass range, as indicated by the line extent (see the text). In 
the right panel, the small symbols connected by solid lines show corresponding results for the Gil simulated galaxy populations. Black 
dashed lines here show the mean luminosity function for all simulated galaxies, rcnormalizcd to fit the satellite data as in the left panel. 
The SDSS results are overplotted as filled triangles for ease of comparison. 



Table 1. Exponent a of the faint-cnd-slope for SDSS luminosity and mass functions 



Range in primary log Ah /Mq 


11.7-11.4 


11.1-11.4 


10.8-11.1 


10.5-10.8 


10.2-10.5 


Luminosity function 
Mass function 


-1.170 
-1.231 


-1.295 
-1.297 


-1.424 
-1.455 


-1.587 
-1.800 


-1.622 
-1.748 



variation of faint-end slope with primary mass seen in the 
SDSS is weak or absent in the simulation data. 

Many previous papers have investigated the 
shape and faint -end slope o f sroup / cluster lum i nosity 



functions (e.g. 

Christlein fc Zabludofj 



Paolillo et ai] 12001 
20031: 



Goto et al.l 

De Propris et al.l 



Andreon. Punzi. fc Gradol 12 005': 'Popcsso et al.ll2005l. 



Zandivarez. Martinez, fc Merchan 2006: Hansen et al.l 



2002 



2002 



2ooe 



200C 



Alshino et al.ll20ld : Ide Filippis et al.ll201ll ). Unfortunatelv 
there is little consensus. Some authors found large differ- 
ences between cluster and field luminosity functions, while 
others found the two to be quite similar. For example, by 
stacking SDSS data around the centres of cluster s detected 
in the Rosat All Sky Survey IPopesso et al.l (|2006l ) obtained 
cluster luminosity functions with a very obvious steepening 
at faint magnitudes. This faint-end upturn appeared to be 
contributed primarily by early-type galaxies. In our most 
massive primary stellar mass bin (11.7 > logjgM* > 11.4), 
there is qualitatively similar behaviour with an upturn at 



Mr -^ — 18 as in lPopesso et al.l (|2006l ) but the steepening is 
much less dramatic in our data than in theirs. In contrast, 
no evidence of an u ptur n at the faint end was found by 
I Alshino etall (|201G| ) and Ide Fihppis et al.1 l|201ll ) in galaxy 
clusters. 



4.2 Stellar mass functions 

Figure |6] is similar to figure [5] but shows stellar mass func- 
tions for satellites both in the SDSS (left panel) and in the 
Gil simulations (right panel with the SDSS data repeated 
as filled triangles). The strong dependence of satellite num- 
ber on primary mass is again evident. The black solid line 
overplotted on the most massive bin is a Schechter func- 
tion fit with characteris tic mass and low-m ass slope fixed 
to the "field" values of iLi fc Wh"it3 |200i). The observed 
satellite stellar mass functions are again steeper than the 
corresponding field function at the low-mass end. The black 
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Figure 6. Similar to figure [5] but showing stellar mass functions for satellites projected within 300 kpc (or 170 kpc) of their primaries, 
both for the SDSS (left panel) and for the Gil simulations (right panel with the SDSS data repeated as filled triangles). As in figurefS] 
the primaries are grouped into five disjoint ranges of log M^^/Mp), as indicate d by the colours. A black solid line in the left panel is a 
renormalized version of the "field" stellar mass function of iLi &: Whitel 1 120091 ) overplottcd on the data for the highest mass primaries, 
while dashed black linos show a power law fit to each mass function estimate. In the right panel the dashed black linos are fits of the 
satellite data to renormalized versions of the stellar mass function of the simulation as a whole. 



dashed lines in the left panel are power-law fits to the obser- 
vational data. The corresponding faint-end-slopes are given 
in the second row of table [l] Here also we have ignored the 
few brightest points in each estimate. As before the extent of 
each dashed black line indicates the points actually used in 
the fit. The steepening of the satellite stellar mass functions 
with decreasing primary mass is even stronger than was the 
case for the satellite luminosity functions. 

In the right panel of figure |6] the Gil points are taken 
from the MS at log M^/M© > 9.5 and from the MS-II at 
lower mass. Fits to the field stellar mass function from the 
simulations are shown as dashed lines and indicate no sig- 
nificant shape difference between the two. Comparison with 
the SDSS data shows that the overprediction of the satellite 
luminosity function for the most massive primaries persists 
in very similar form in their stellar mass function. The un- 
derprediction found for the two lower primary mass ranges 
has gone away, however, indicating that the discrepancy 
was due primarily to the colours of the simulated satellites, 
rather than to their stellar masses. The simulation does not 
reproduce the clear steepening of the SDSS satellite mass 
functions with decreasing primary mass. Any effect in this 
direction is very weak. 

To investigate further the origin of the observed steep- 
ening with decreasing primary mass, we take all spectro- 
scopic galaxies brighter than 16.6 in the r-band, without 
applying any isolation criterion, and calculate the mass func- 
tion of surrounding satellites within 300 (or 170) kpc in ex- 
actly the same way as for our isolated primaries. Results 
are shown as solid curves in figure IX6l In this case there is 



no steepening with decreasing primary mass - the different 
colour curves are more or less parallel to each other and to 
the field stellar mass function. Apparently, the steepening 
is caused by our isolation criteria - isolated galaxies have 
fewer lower mass neighbours than typical galaxies. The sup- 
pression is stronger for more massive neighbours, steepen- 
ing the satellite mass functions of isolated primaries. Notice 
that for primaries with logM*/MQ > 10.8, the mean num- 
bers of low-mass companions are similar for isolated and 
for typical objects, but that lower mass primaries have sub- 
stantially fewer such companions if they are isolated. This 
is because many of the lower mass non-isolated "primaries" 
are, in fact, satellites themselves, and their low-mass "com- 
panions" are fellow satellites within the larger system. It 
is quite interesting that the abundance of relatively small 
satellites is strongly affected by the presence or absence of a 
nearby galaxy comparable in luminosity to the primary, par- 
ticularly since the observed trends are not fully reproduced 
by the simulation. This suggests that dwarf galaxy forma- 
tion may be influenced by nearby giants in a way which the 
simulation does not represent. 

It is interesting to see whether satellite galaxy popu- 
lations depend on the colour of the primary galaxy as well 
as on its stellar mass. Figures [7] and |8] show the satellite 
mass functions surrounding red and blue primaries respec- 
tively, where the primary populations have been split at the 
colours indicated in figures [l] and [2] The black dashed lines 
in the left panels are power-law fits, but in each case the 
slope is fixed to be that found for all primaries in the rele- 
vant stellar mass bin (table [l}. In the right panel of figure [8] 
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Figure 7. Similar to figure [6] but for red primaries. See figures [T] and [2] for tlie colour cuts separating red and blue primaries. 
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Figure 8. Similar to figure [B] but for blue primaries. In the right panel, the red dots connected by a curve stop at log^gM* 
because the Millennium-II Simulation has fewer than eight blue primaries more massive than log M-^/Mq = 11.4. 



9.5, 



the red curve stops at logAf*/MQ — 9.5 because the MS-II 
contains fewer than eight blue primaries more massive than 
log Mi,/MQ = 11.4 and the MS population is affected by nu- 
merical resolution at lower satellite mass. For red primaries 
the SDSS and Gil data agree quite well, apart from the 
slope discrepancy and a residual overprediction of the abun- 
dance of faint satellites around high-mass primaries. For 



blue primaries, the differences are bigger but in the high- 
est mass bin this could reflect the small number of SDSS 
primaries and the correspondingly large observational error 
bars. The discrepancy for primaries in the stellar mass range 
10.1 > log M^/Mq > 10.8 is smaller but more significant, 
particularly since there is good agreement in this mass range 
for red primaries. Interestingly, a comparison of figures [7] 
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and [8] shows that the amphtude of the satelhte luminosity 
function is higher around red primaries than around blue 
primaries of the same stellar mass. We analyze this result in 
more detail in the following subsections. 



4.3 Satellite abundance as a function of primary 
stellar mass and colour 

In order to better display how the abundance of satellites 
depends on the stellar mass and colour of the primary, we use 
the power-law fits shown as dashed lines in figures[Sl [7]and[S] 
to predict the mean number of satellites per primary in the 
stellar mass range 10.0 > log Mi^/Mq > 9.0 and at projected 
separation Vp < 300 kpc (rp < 170 kpc for the lowest mass 
primaries). This has the advantage of producing a robust 
measure of satellite abundance which is little affected either 
by selection-induced cut-offs (most important for low-mass 
and red primaries) or by incompleteness (most important 
for high-mass and blue primaries). The results are shown 
in the left panels of figure [^ where black dots and lines 
give results for all primaries, while red and blue dots and 
lines give results for red and blue primaries, respectively. 
The top panel presents results for the SDSS and the lower 
panel results for the Gil simulations. The SDSS result for 
all primaries is repeated in the lower panel, showing that the 
simulation overpredicts the number of satellites in this mass 
and projected radius range both for the highest mass and 
for the lowest mass primaries. Note that because the low- 
mass slopes differ in simulation and observation, the result 
for low-mass primaries depends on the satellite mass range 
chosen for the comparison. 

At high mass the black and red curves in figure |9] are 
close to each other, reflecting the fact that the fraction of red 
primaries is large (see figure[31). At the highest mass, the blue 
curve indicates a consistent number of satellites around blue 
primaries, although with considerable uncertainty because 
such primaries are rare. At somewhat lower mass, however, 
blue primaries have significantly fewer satellites than red 
primaries of the same stellar mass, both in the SDSS data 
and in the simulation. This is a primary result of our paper. 
The effect is a factor of two to three in satellite abundance 
for primaries with log M^/Mq ~ 11. In the SDSS data there 
is some indication that the the colour dependence may get 
smaller again for lower mass primaries, but this does not 
happen in the simulation, where there is still more than a 
factor of two difference for A'U ~ 2 x 10^° Mq. Overall, the 
differences appear somewhat larger in the model than in the 
real data. 

The cause of this effect in the Gil simulation is easy 
to track down. In the right-hand panels of figure [9] we plot 
histograms of host halo mass for isolated galaxies as a func- 
tion of their stellar mass and colour. (As shown in figure [31 
almost all isolated galaxies in the simulation are the central 
galaxies of their haloes.) For all except the highest stellar 
mass range, red primaries have significantly more massive 
dark haloes than blue ones. The shift between the peaks of 
the two distribution is an order of magnitude for primaries 
with 11.4 > log Afj,/M0 ~ 11.1, dropping to a factor of two 
for 10.5 > logMt/Af© ^ 10.2. In the simulation red pri- 
maries have more satellites because they live in more mas- 
sive haloes. A direct indication that the same may hold for 
real galaxies comes from the galaxy-galaxy lensing study of 
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Figure 10. Satellite mass functions within rp = 300 kpc, split ac- 
cording to halo mass, for dark haloes in the Gil galaxy formation 
simulation. The colours refer to different ranges of logMjiaio/Af© 
as indicted in the legend. For each halo mass range, satellite stel- 
lar mass functions are shown separately for haloes with red (solid 
lines) and blue (dashed lines) central galaxies. 



iMandelbaum et al.1 (120061 1. By combining their SDSS lens- 
ing data with HOD modeling, these authors concluded that 
red galaxies have more massive haloes than blue ones for 
logAf*/M0 > 11. At lower central galaxy masses their re- 
sults appear consistent with no offset, although the error 
bars are large (see their figure 4). 

The fact that both the number of satellites and the mass 
of the associated dark halo depend not only on the stellar 
mass of the primary galaxy but also on its colour contra- 
dicts the assumptions underlying many HOD or abundance 
matching schemes for interpreting large-scale galaxy clus- 
tering. Such a dependence could be included in more com- 
plex versions of at least the former, but would require addi- 
tional parameters and additional ob s ervational data to con- 
strain them (e.g . Simon et al.l l2009l : IRoss 



2009H . 



Brunneij[2009l : 



ISkibba fc Shetll l2 

Within the simulation it is possible to check whether 
halo mass is the only factor responsible for the difference in 
satellite abundance between red and blue primaries. In fig- 
ure [TD] we present stellar mass functions for satellites of iso- 
lated galaxies as a function of the halo mass of the primary. 
As before these are compiled for satellites projected within 
rp = 300 kpc. The halo mass ranges for this plot are chosen 
to correspond roughly to the primary stellar mass ranges in 
previous figures. For each halo mass range, the mass func- 
tions are also split according to the colour of the primary 
galaxy, with solid and dashed lines referring to results for red 
and blue primaries respectively. At low satellite mass there 
is excellent agreement between the solid and dashed curves, 
indicating that satellite abundance does not depend on cen- 
tral galaxy colour at fixed halo mass. For the most massive 
haloes there are few blue primaries in the MS-II, so the red 
dashed curve is quite noisy below logM^^/AlQ — 9.5. For 
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Figure 9. Left: Mean number of satellites in the stellar mass range 10.0 > log M^/Mq > 9.0 as a function of primary stellar mass. 
Blaek, red and blue points refer to all, to red and to blue primaries, respectively. For the points at lowest primary mass, satellite counts 
were accumulated within 170 kpc, whereas for all other primary masses they were accumulated within 300 kpc. The top panel gives 
observational results for the SDSS while the bottom panel gives corresponding results for the Gil galaxy formation simulations. The 
SDSS result for all primaries is re-plotted as a dashed curve in the bottom panel in order to facilitate comparison. Right: Host halo mass 
distributions for red and blue primaries in the simulated catalogues, split into the same primary stellar mass ranges as in the left-hand 
plots. 



massive satellites there are obvious discrepancies between 
the dashed and solid curves, but these result from our sam- 
ple definition. At given halo mass, red primaries have smaller 
stellar masses and so substantially lower luminosities than 
blue ones. Our isolation criteria then imply a correspond- 
ingly lower upper limit on the stellar mass of satellites for 
the red primaries. The effect is largest for the lowest mass 
haloes. 

I Sales et al. 



||2007^ used the galaxy formation simu- 



lation of ICroton et al.l (|2006l '). also based on the Millen- 
nium Simulation, to study the relation between primary 
luminosity/stellar mass and satellite velocity dispersion, 
which should be a good diagnostic of halo mass. Although 
they found a strong dependence of velocity dispersion on 
galaxy colour at fixed primary luminosity, this dependence 
almost vanished at fixed primary stellar mass. This ap- 
pears to contradict our results from the Gil simulation. 
At log M^/M0 ~ 11, where we find the biggest difference 
in satellite abundance between red and blue primaries, al- 
most a factor of three, the difference in velocity dispersion 
in their figure 13 is at most 20%, corresponding to a factor 
of at most 1.7 in halo mass (since Mh oc a^). The discrep- 
ancy could result from the different isolation criteria adopted 
in the two studies, from departures from a straightforward 
relation between 3-D velocity dispersion within r2oo (the 
quantity considered bv I Sales et al.l (|2007|)), halo mass and 
projected satellite count within 300 kpc (the quantity con- 
sidered here), or from differences between the two galaxy 
formation models. 

So far we have characterized the abundance of satellites 
by the count within Vp = 300 kpc (or r^ = 170 kpc for the 
lowest mass bin in figure ^ In order to better understand 



the relation with halo mass, it is useful instead to consider 
the count within the virial radius of the haloes, which we 
define as r2oo , the radius of a sphere within which the mean 
mass density is 200 times the critical value. We obtain such 
a measure using the equation, 

M{r < r2oo) 



N{r < r-2oo) = Np{Tp < 300 kpc) x 



(4) 



Mp{rp < 300 kpc) 

where Mp and Np are projected quantities and M and N 
the corresponding 3-dimensional quantities. The projected 
abundance Np{rp < 300 kpc) is taken directly from fig- 
ure [9] and we calculate the mass ratio on the rhs of equa- 
tion |4] assuming an NFW profile . The halo radius r2oo is 
known for each simulated galaxy, but can only be inferred 
indirectly for the SDSS objects. We assume that the mean 
halo mass at given stellar mass is the same in the observa- 
tions as in the Gil simulations. The mass ratio also depends 
w eakly on hal o conc entration which we take from the model 
of lZhao et al] (|2009l ). 

For each primary stellar mass bin in figure [5] we can 
then estimate a mean halo mass (M2oo)- We take the mean 
satellite number per primary given by equation |4] and di- 
vide by this mean halo mass to obtain the abundance 
per unit halo mass of satellites in the stellar mass range 
10.0 > log M*/M0 > 9.0. Finally, we divide this quantity 
by the abundance per unit (total) mass of galaxies in this 
same stellar mass range in the Universe as a whole, taken 
from iLi fc Whitel (|2009l ) for the SDSS data and from the 
Gil simulation as a whole for the mock data. The result is 
a measure of the formation efficiency of low-mass galaxies 
as a function of their present environment, as characterized 
by mass of the halo in which they live. Figure [TT] shows this 
efficiency as a function of mean halo mass for all primaries 
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Figure 11. Mean number of satellites in the stellar mass range 
10.0 > logM*/MQ > 9.0 per unit total halo mass, relative to 
the mean abundance per unit total mass of such galaxies in the 
Universe as a whole. Isolated primary galaxies arc grouped into 
the same five stellar mass bins used in figure |9] Results for the 
Gil simulations and for the SDSS are shown by solid and dashed 
curves respectively. The mean halo mass has been calculated di- 
rectly for each bin in the simulation data. Each SDSS primary is 
assigned a halo mass using the simulation relation between mean 
halo mass and primary stellar mass for all primaries. This plot 
hence shows the efficiency of low-mass galaxy formation as a func- 
tion of present-day halo mass in units of the overall efficiency in 
the Universe as a whole. 



in the five stellar mass bins of figure [9] There are three sig- 
nificant points to take from this plot: (i) the formation effi- 
ciency of low-mass satellite galaxies varies rather little with 
the mass of the halo in which the galaxies are found today; 
(ii) in massive haloes this efficiency is about 50% larger than 
the efficiency for forming such galaxies in the universe as a 
whole (remember that, globally, about half the galaxies in 
this stellar mass range are satellites); (iii) finally there is fair 
agreement between the formation efficiencies in the simula- 
tion and in the real universe, although this is in part due 
to our use of the simulation to assign halo masses to the 
observed galaxies. 



5 SATELLITE COLOUR DISTRIBUTIONS 

So far we have studied the abundance of satellites as a func- 
tion of the stellar mass and colour of their primary and as 
a function of their own luminosity and stellar mass. In this 
section we study how the colours of satellite galaxies depend 
on the properties of their primaries. Figures [T^ and [T51 show 
cumulative colour distributions for satellites in the SDSS 
and in the Gil catalogues respectively, as a function of the 
stellar mass and colour of their primary. The distributions 
are for satellites in two different stellar mass ranges, as in- 
dicated by the labels above the relevant panels, and refer 
to all satellites projected within 300 kpc. The top, middle 
and bottom panels in each column refer to satellites of all, 
of red and of blue primaries respectively, while the differ- 
ent colours of the curves in each panel encode the stellar 
mass range of the primary galaxies. The black curves which 
repeat in all the panels of each column give the colour dis- 
tributions for field galaxies in the same stellar mass range 
as the satellites (calculated for SDSS from all galaxies in the 



NYU-VAGC with r < 17.6, and for Gil from all galaxies 
within the simulation volume). The dashed horizontal line is 
merely a reference to facilitate identification of the median 
colour. 

A number of systematic trends are evident in these 
plots. Concentrating first on the observational results in 
figure 1121 we see that more massive primaries have redder 
satellites within 300 kpc (in every panel the curves are or- 
dered cyan- blue-green-red from top to bottom), that low- 
mass satellites are bluer than high-mass ones (the curves 
in the right panels are always bluer than the correspond- 
ing curves in the left panels), satellites are systematically 
redder than field galaxies of the same mass, except possibly 
for the lowest mass primaries (the coloured curves almost 
always lie below the corresponding black curves), and red 
primaries have redder satellites than blue primaries of the 
same stellar mass (every coloured curve in the lowest panels 
is bluer than the corresponding curve in the middle panel). 
This l ast trend is the "galact ic conformity" effect pointed 
out bv lWeinmann et all (|2006l '). 

If we now compare with the simulation results in fig- 
ure [13] we see that the same four systematic trends are 
present. More massive primaries have redder satellites; lower 
mass satellites are bluer; satellite galaxies are redder than 
field galaxies of the same stellar mass; and red primaries 
have redder satellites than blue primaries of the same stel- 
lar mass. However, there is an obvious discrepancy in that 
simulated satellites are systematically redder than observed 
satellites. This is true for all primary and satellite masses, 
but is particularly marked for lower mass and red primaries, 
and for lower mass satellites. Clearly, the theoretical model 
of Gil suppresses star formation much more effectively in 
such satellites than is t he case in the real uni verse. This 
echoes the concl usions of Weinmann et al.l ( 20061 ') about the 
earlier models of lDe Lucia fc Blaizoj (20071). The excessive 
reddening of the simulated satellite population reduces but 
does not eliminate all the other trends mentioned above. 

The galactic conformity phenomenon has been 
discus sed in a number o f previous publi c ations 
fe.g. IWein mann c t all l2006l: lAnn. Park, fc Choil 12008 ' 



Kauffmann. Li. fc Heckmanl |2010| : IPrescott et ahl 120111 ). 



Weinmann et al.l ( 20061 ) showed that, among groups of 



given luminosity (which they considered a proxy for halo 
mass), those with an early-type central galaxy have a 
larger fraction of early-type satellites. They considered 
several physical processes which might be responsible for 
this (halo and/or galaxy mergers, ram-pressure stripping, 
strangulation, harassment...), focusing on whether these 
processes could alter galax y morphology. Howeve r, as 
discussed in some detail by iKauffmann et al.l (|2004 ) and 
re-emphasized in the cont ex t of galactic conformity by 
IKauffmann. Li. fc Heckmanl (|2010|), it is important to 
separate star formation activity, stellar mass and galaxy 
structure when analyzing the infiuence of environment on 
galaxy properties. Typical classifications into "early" and 
"late" types mix aspects of all of these. The conformity 
effects we see here are for given central galaxy stellar mass 
(rather than luminosity), are within a fixed projected radius 
(300 kpc), and refer specifically to the colours of satellites 
and primaries. It seems possible that they could be due at 
least in part to the tendency for red centrals to have more 
massive haloes than blue ones, together with the trends 
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for satellites to get redder with increasing halo mass and 
decreasing r/r2oo- 

This can be checked directly for the simulated galaxy 
catalogues of Gil. In figure [14] we again show cumulative 
satellite colour distributions for two different satellite mass 
ranges and for all, for red and for blue primaries, but now in 
bins of host halo mass rather than of primary stellar mass. In 
this plot, the colour distributions are calculated for all satel- 
lites projected within the halo virial radius rather than for 
a fixed projected radius of 300 kpc. The halo mass ranges 
have been chosen to correspond approximately to the pri- 
mary stellar mass ranges we have been using in previous 
plots. For the two high-mass bins, the colour distributions 
show no significant dependence on the colour of the central 
galaxy, but for the two low-mass bins the dependence on pri- 
mary colour, while smaller than in figure [T3l is clearly still 
present. Thus there must be physical processes in the model 
which contribute to the "galactic conformity" phenomenon 
in addition to those which result in red primaries having 
more massive haloes than blue primaries of the same stellar 
mass. 

We have analyzed our simulation to identify candidates 
for these additional processes. In figures [TS] and [TS] we plot 
distributions for the two which show the strongest trends. 
Figure [T5] shows the distributions of black hole mass and 
of hot gas mass for haloes in these same four mass ranges, 
split according to the colour of the central galaxy. Figure [T6l 
shows similar plots for the cumulative distributions of in- 
fall redshift (defined as the redshift when a satellite last 
entered the virial radius r2oo of the main progenitor of its 
halo) for satellites with stellar mass above W^Mq. In all 
cases the distributions depend strongly on central galaxy 
colour for haloes in the two lower mass ranges where galac- 
tic conformity effects are substantial, but at most weakly 
for higher mass haloes where these effects are absent. For 
logA/h/M0 < 12.8, haloes with red central galaxies have 
more hot gas, more massive central black holes and ear- 
lier satellite infall redshifts than haloes of the same mass 
with blue central galaxies. Notice that the black hole mass 
distribution for blue centrals is bimodal in the lower left 
panel of figure 1151 Clearly, there is a transition in this halo 
mass range between blue centrals with high-mass black holes 
(about 0.3% of their stellar mass, as in the upper panels) 
and blue centrals with low-mass black holes (roughly an or- 
der of magnitude smaller as a fraction of stellar mass, as 
in the lower right panel). The hot gas fractions are similar 
to the overall cosmic baryon fraction (~ 17%) for haloes 
with red central galaxies and about a factor of two smaller 
in the lower mass haloes with blue centrals. The transition 
halo mass is approximately the mass at which both obser- 
vational and simulated samples of isolated galaxies become 
dominated by red objects (see figure [3]). 

Given the modeling assumptions of Gil, these system- 
atic differences mean that, for log M^/Mq < 12.8, haloes 
with red central galaxies have gas cooling rates which are 
twice but "radio mode" heating rates which are 20 times 
those of similar mass haloes with blue central galaxies. Thus, 
their central galaxies are clearly red because feedback has 
quenched their growth, and this is why their stellar masses 
are smaller than those of the blue centrals which have con- 
tinued to grow to the present day. The redness of the satel- 
lite population in haloes with red centrals can be traced to 
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Figure 12. Cumulative colour distributions for satellites pro- 
jected within 300 kpc of their primary as a function of primary 
stellar mass (indicated by line colour) for two ranges of satellite 
mass (left and right columns with the range indicated by the la- 
bel above each column) and for all, for red and for blue primaries 
(upper middle and lower panels in each column) . Black lines show 
the cumulative colour distribution for a volume-limited sample of 
field galaxies derived from the SDSS DR7 spectroscopic sample 
with r < 17.6. 



the facts that the satellites are accreted earlier and orbit 
through a denser hot gas medium. Thus both tidal and ram 
pressure stripping processes are more effective, and the satel- 
lites have had longer to exhaust any remaining star-forming 
gas. Finally, the earlier assembly indicated by the higher 
satellite infall redshifts presumably explains why these par- 
ticular haloes have red centrals, bigger black holes and more 
hot gas. 

Thus, at least in the models, it is clear why there is 
galactic conformity at given halo mass. The more easily ob- 
servable galactic conformity at given central stellar mass 
is predicted to be stronger as a result of the combination 
of these effects with the tendencies for for red centrals to 
have more massive haloes at given stellar mass than blue 
ones, and for more massive haloes to have redder satellites. 
Yet stronger conformity effects are expected at given central 
galaxy luminosity, since blue central galaxies have smaller 
stellar masses (and thus even smaller halo masses) than red 
central galaxies of the same luminosity. 
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Figure 13. Similar to figure [T2l but for the simulated galaxy cat- 
alogues of Gil. The black lines here are the colour distributions 
for all galaxies in the specified stellar mass ranges within the full 
simulation volume. 
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Figure 14. Colour distribution of satellites projected within the 
halo virial radius r200 of all primaries, of red primaries and of 
blue primaries (top, middle and bottom rows) as a function of 
host halo mass (indicated by line colours corresponding to the 
ranges of log Mi^/Mq given in the legend) in the simulated galaxy 
catalogues of Gil. Results are again shown for two ranges of 
satellite stellar mass. 



6 SUMMARY AND CONCLUSIONS 

We have used a photometric catalogue of SDSS/DR8 galax- 
ies brighter than r = 21 to study the satellite populations 
of 41271 isolated galaxies with r < 16.6 selected from the 
SDSS/DR7 spectroscopic catalogue. In particular, we have 
studied how the abundance of satellites as a function of 



luminosity, stellar mass and colour depends on the stellar 
mass and colour of the central galaxy. Our study differs 
from other recent SPSS-based studies of satellite g alaxies 
(jLares. Lambas. fc DominguejuOllI : IGuo et al.ll201ll ) in the 
size of the sample analyzed, in our primary focus on sys- 
tematics as a function of stellar mass, and in our detailed 
comparison with the predictions of simulations of the evo- 
lution of the galaxy population in the concordance ACDM 
cosmology. In general, our results confirm and extend those 
obtained earlier, but the comparison with simulations allows 
us to identify the likely physical cause of most of the effects 
we see, and to isolate those which do not have a natural 
explanation within our current theory of galaxy formation. 
Our observational samples and analysis procedures al- 
low us to measure the properties of the satellite popu- 
lation in an unbiased way down to absolute magnitudes 
Mo.i^ ~ — 14 and stellar masses logM^/MQ ~ 8. Our main 
observational conclusions are as follows; 

• Satellite luminosity and stellar mass functions have 
shapes consistent with those of the general field galaxy 
population only around the highest stellar mass primaries 
logAf*/Af0 > 11.4. These are all brightest cluster galaxies. 
For lower mass primaries, these functions become progres- 
sively steeper, even after accounting for the bright-end cut- 
off induced by our isolation criteria. This steepening is more 
marked for the stellar mass functions than for the luminos- 
ity functions because observed satellites get bluer as their 
stellar mass decreases. 

• The mean abundance of satellites increases strongly 
with primary stellar mass, approximately as expected if the 
number of satellites is proportional to dark halo mass. 

• For logM*/MQ > 10.8, red primaries have more satel- 
lites than blue primaries of the same stellar mass. The ef- 
fect exceeds a factor of two for log M^t/Mp ~ 11.2. T his is 
reminiscent of the result of iMandelbaum et al.l (|2006l ) who 
showed that at high stellar mass, red central galaxies have 
more massive haloes than blue ones. This trend could in part 
be due to colour dependent errors in deriving stellar masses 
from the photometry, but such errors would need to be quite 
large and to depend on primary mass. 

• Satellite galaxies are systematically redder than field 
gala^xies of the same stellar mass except around blue pri- 
maries with log Mi,/MQ < 10.8 where the satellites can have 
similar colours or even be systematically bluer than the field 
(i.e. the galaxy population within a large representative vol- 
ume). 

• The satellite population is systematically redder around 
more massive primaries, for more massive satellites and 
around red primaries. The first effect reflects the fact that 
cluster galaxies are systematically redder than fleld galax- 
ies, the second echoes the trend found in the general field, 
and the third is the gala ctic conformity effect pointed out by 
IWeinmann et al.l (|2006l ) but measured here for fixed central 
stellar mass rather than fixed central luminosity. 

We used criteria directly analogous to those employed 
on the SDSS to construct an isolated galaxy sample from 
the 2 = outp ut of the publicly available galaxy formation 
simulations of IGuo et al.l (|201ll . Gil). These are based on 
the Millennium and Millennium-II Simulations. The mock 
catalogue contains similar magnitude, stellar mass, colour 
and position/velocity information to the real catalogue, but 
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also contains information about dark haloes and the loca- 
tion of the galaxies within them. Based on the mock sample, 
we conclude that ~ 98% of our isolated galaxies are the cen- 
tral objects of their dark haloes. Both in the SDSS and in 
the mock catalogue, the distributions of intrinsic proper- 
ties for the isolated and parent populations are very similar. 
Only the colour distributions shift slightly, with the isolated 
galaxies being systematically bluer than the full population. 
A detailed comparison of the mock and real samples leads 
to the following conclusions 

• At all primary masses, the luminosity and stellar mass 
functions of Gil satellites are quite similar both in shape 
and in normalization to those measured for SDSS. However, 
in the simulation there is only a weak tendency for the satel- 
lite functions to be steeper than those of the field, or to be 
steeper for lower mass primaries. This disagrees with the 
SDSS where the steepening with decreasing primary mass is 
quite marked. 

• In the mock catalogues the abundance of satellites in- 
creases with primary stellar mass almost in proportion to 
mean halo mass. For high-mass haloes, the abundance per 
unit mass of satellite galaxies is about 1.5 times the value 
for the universe as a whole. 

• For log M*/M0 < 11.4, red simulated primaries have 
more satellites then blue ones of the same stellar mass. The 
effect is similar in strength to that seen in SDSS but con- 
tinues to lower primary mass. In the simulation it is due 
entirely to red primaries having more massive haloes than 
blue ones. At fixed halo mass the abundance of faint satel- 
lites is independent of the colour of the primary, but red 
primaries have lower stellar masses because of their trun- 
cated star formation histories. 

• Satellite galaxies in the Gil simulations are systemat- 
ically redder than in the SDSS. The effect is particularly 
marked for lower mass satellites and around lower mass pri- 
maries. The modeling improvements introduced by Gil to 
address this issue are apparently insufHcient to fully solve 
it. Star formation is still terminated too early when galaxies 
become satellites. 

• Despite this overall shift, the colours of simulated satel- 
lites depend on the colour and stellar mass of their primary 
and on their own stellar mass in very similar ways as in 
the SDSS. Satellites are systematically redder if they are 
more massive, if their primary is more massive, and if their 
primary is red. The first trend echoes that found for field 
galaxies, while the second reflects the fact that more massive 
haloes contain redder satellites. The third trend, the galactic 
conformity effect, is caused by redder primaries having more 
massive haloes at fixed stellar mass, and denser hot gas at- 
mospheres (hence more effective ram-pressure stripping) at 
fixed halo mass. 

Satellites go red in the Gil simulations because they 
lose their source of new cold gas. Once they fall within 
the virial radius of their host, their hot gas reservoirs are 
gradually removed by tidal and ram-pressure stripping and 
no new material is added by infall. Star formation uses up 
the remaining cold gas a nd then switches off. C learly this 
happens too quickly (see IWeinmann et al.l (|201ll ) and also 
IWang et al.l ((2007]) for an explicit demonstration of how 
lengthening the relevant timescales can cure the problem) 
so improving the model will require changing the star for- 




Figure 16. Infall time distributions for satellites in the same four 
ranges of halo mass used in figure [14] and again split according to 
the colour of the central galaxy. 



mation assumptions to increase the time to gas exhaustion 
(for example, by removing the threshold gas surface density 
for star formation) or providing new sources of fuel (for ex- 
ample, by including the gas return from stellar evolution). 
Once this is fixed, it seems likely that the other trends of 
satellite colour with environment will be well matched. 

The other clear discrepancy between the SDSS data 
and the simulation is the steepening of satellite mass and 
luminosity functions as one goes to fainter primaries. This 
is a relatively strong e ffect in the real data (and is visible 
also in the analysis of IGuo et al.l (|201H )) but does not oc- 
cur at a significant level in the simulations. Thus, it must 
refiect not merely the statistics of hierarchical clustering in 
a ACDM cosmology, but in addition some difference in the 
star-formation histories of satellite galaxies living in differ- 
ent mass host haloes. This may be related to the discrepancy 
noted by Gil - comparisons with high-redshift data suggest 
that low-mass galaxies form too early in their simulations - 
but this can only be checked by further simulation work. 
It is clear that the detailed and relatively precise statisti- 
cal information provided by large-sample studies of satellite 
galaxies is useful for testing and refining our understanding 
of how galaxies form. In a follow-up paper we will extend 
our current study by considering satellite galaxy properties 
as a function of distance from the primary. 
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Figure 15. Left: Central black hole mass distributions for the same four ranges of halo mass used in figure [141 and split according to the 
colour of the central galaxy. Right: Hot gas mass distributions for the same four sets of haloes and again split according to the colour 
of the central galaxy. 
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APPENDIX A: DATA SYSTEMATICS AND 
COMPLETENESS 

Al SDSS Photometry Systematics 

As mentioned by iMandelbaum et al.l (120051. 2006 ) and also 
discussed in the SDSS/DR8 paper JAihara et al.ll201ll '). the 
magnitudes of faint objects close to brighter galaxies are 
affected by systematic problems related to sky subtraction 
and deblending. This results in an apparent over- or under- 
abundance of faint galaxies close to bright ones. This is po- 
tentially a critical problem for our study, since we evaluate 
the mean number of satellites surrounding some set of pri- 
maries by counting all faint neighbours in the photometric 
catalogue and then subtracting the number of background 
objects that should be projected onto these same regions "at 
random" . We therefore test explicitly for such affects in the 
DR8 photometry we will use. 

We take our set of isolated primary galaxies and use 
the techniques described in the main body of our paper to 
count "red" neighbours brighter than r = 21 in a set of 
annuli of fixed angular scale. We also make similar counts 
using a fake photometric catalogue built by randomizing the 
positions of the actual photometric galaxies within the sur- 
vey footprint. We then accumulate both sets of counts for 
each annulus and for all primaries within each of the five 
stellar mass bins used in figure 1, and we divide the real 
counts by the random counts to get a normalized profile. 
Error bars are estimated by bootstrapping the sample of 
primaries. Our definition of "red" is that the colours of the 



photometric galaxy should be too red for it to be at the 
same redshift as its primary. We estimate a stellar mass M-,, 
for the photometric galaxy assuming it to be at the redshift 
of the primary (see section |3}. We then consider any photo- 
metric galaxy with "'^(g - r) > 0.0321og^(,M* -f 0.73 to be a 
background object. (This is based on the upper envelope of 
a rest-frame colour versus stellar mass scatter plot for SDSS 
galaxies with measured redshift.) Our normalized profile is 
thus expected to be unity to within its statistical errors. The 
actual results are shown in the left panel of figure I Ali tor our 
five primary stellar mass bins. 

All the profiles go nicely to unity at angular separa- 
tions of an arcminute or more. This shows that our colour 
cuts have indeed isolated a population of background galax- 
ies which has no significant spatial correlation with the pri- 
maries. On the other hand, there are substantial deviations 
from unity at smaller scale which can have either sign and 
are similar for the different primary mass ranges. These most 
likely refiect photometric errors caused by the effect of the 
bright galaxy on the local sky background estimate used 
when meas uring the apparent magn itude of its faint "red" 
neighbour (JMandelbaum et al]|2006h . At r = 21 the loga- 
rithmic slope of the integral r-band number counts is 0.34 
dex per magnitude, so a change in the apparent surface den- 
sity of background objects by a factor of 1.5 requires an ef- 
fective shift of about 0.5 mag in the apparent magnitude 
scale[j Note that the typical angular size of the primaries 
(characterized by (690), where ^90 is the angular scale con- 
taining 90% of the Petrosian luminosity, and indicated by 
the arrows in the individual panels of figure IXTj) is quite sim- 
ilar for each of our stellar mass bins and is comparable to 
the separations where we see substantial systematic effects. 

These problems will affect the SDSS photometry both 
of true satellites and of background galaxies, presumably by 
similar amounts. We will neglect the shifts in satellite mag- 
nitudes in the following since they are only substantial at 
separations below 20 arcsec, and so will affect very few ob- 
jects. A more serious issue is whether the apparent excess 
(or deficit) of background galaxies will cause our background 
subtraction to fail, since this assumes the number density of 
background galaxies to be statistically uniform and inde- 
pendent of position on the sky relative to the primary. We 
estimate the size of the problem as follows. We construct a 
random sample by taking the full photometric catalogue and 
randomizing the positions of the galaxies within the survey 
footprint, keeping all their other properties fixed. We then 
accumulate the total count A^tot of random galaxies with 
r < 21 in each radial annulus for all the primaries in each 
mass bin. We multiply A^tot by the factor A*'rcd/A*'rcd,rand — 1. 
(This is the quantity on the y-axis in the left panel of fig- 
ure I All minus unity.) The result is our estimate A'^prob of the 
number of apparent satellites brighter than r = 21 result- 
ing from incorrect background subtraction. This systematic 
error will cause us difficulties if Aprob is comparable to the 
total number of satellites. 

In the right panel of figure [XT] we show the ratio of this 
estimated number of spurious satellites to the total number 



^ Gravitational lensing could also change the apparent density of 
background objects, but simple estimates suggest the eflect is too 
small to be relevant. 
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Figure Al. Left: The ratio of the number profile of "red" background galaxies with r < 21 around our primary sample to a similar profile 
for the same primaries made after randomizing the positions of the photometrie galaxies within the survey footprint. The five panels refer 
to our five ranges of logM*/MQ as indicated by the labels. The arrows show (690), the mean angular size of the primaries in each panel. 
Right: The ratio of our estimate of the number of "false" satellites brighter than r = 21 resulting from systematic errors in background 
subtraction to our estimate of the total number of satellites above this same apparent magnitude limit. The error bars are obtained by 
propagation of those in the left panel, together with bootstrapping the primary sample when estimating the satellite+background counts. 



of satellites we find with r < 21 when the techniques dis- 
cussed in section [3] are applied to the real data. It seems 
that this systematic error is not a problem for the four mas- 
sive primary bins where jAfprobj/A'sat is less than 20% on all 
scales. In the least massive bin, however, |A'prob|/A^Bat seems 
to be around 40% to 50% at Vp > 30kpc, though with pretty 
large error-bars. For this primary mass range the number of 
satellites is very small, and thus even a small misestima- 
tion of the number of background galaxies can cause a rel- 
atively big effect. Looking at the scale range corresponding 
to Vp > 30kpc in the left panel of that bin (this is roughly 
6 > 30" for the relevant mean redshift z ~ 0.06) we find 
that the background density modulation required is smaller 
than 10%, corresponding to an apparent magnitude error of 
~ 0.1 mag. A shift of this order in magnitude due to back- 
ground subtaction problems does not seem implausible. 



We conclude that our results are insensitive to these 
small-scale photometric problems except for the lowest stel- 
lar mass range we consider, where satellite counts may be 
significantly affected. We do not attempt to correct this sys- 
tematic because we are unable to estimate sufficiently ac- 
curately how the effect depends on the photometric prop- 
erties (apparent magnitude, colour, profile shape...) of the 
brighter and fainter images. One might hope that the effects 
would become smaller at brighter apparent magnitudes, but 
we have been unable to demonstrate this convincingly. We 
note that we made similar tests using DR7 photometry, find- 
ing substantially larger effects than presented here for DR8. 
This is our primary reason for choosing DR8 photometry 
when searching for satellites, even though our primary sam- 
ple was defined using DR7 data. 
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Figure A2. Scatter plot comparing SDSS and COSMOS i-band 
apparent magnitudes (including both galaxies and stars). Panels 
on the top and on the right show the fraction of COSMOS objects 
unmatched in SDSS and of SDSS objects unmatched in COSMOS, 
respectively. 



A2 Completeness in SDSS 

In this paper we study the properties and the distribution of 
satellite galaxies down to faint apparent magnitudes where 
identification, classification and magnitude calibration could 
be affected by a variety of observational factors such as see- 
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ing, sky brightness and extinction, in addition to the issues 
of photometric accuracy in the neighborhood of brighter ob- 
jects which we addressed in the previous section. In this 
section we focus on estimating quantitatively tire fraction 
of galaxies that might have been missed in the SDSS cat- 
alogues. We test the accuracy of star-galaxy separation in 
section IA31 of our background galaxy subtraction methods 
in IA4I and the influence of additional observing factors in 
section IA5I 

We investigate the completeness of the SDSS/DR8 pho- 
tometric catalogue using the much deeper HST data avail- 
able for the 2 squ.deg. COSMOS survey. We cross-match 
galaxies and stars in the COSMOS field [f| to galaxies and 
stars in SDSS/DR8 and examine the failed matches in both 
directions. The COSMOS survey was initially conducted 
by HST in broad 7-band using the Advanced Camera for 
Surveys (ACS) Wide-Field Channel (WFC) detector. The 
current COSMOS database also includes later observations 
from a variety of telescopes. It reaches a limiting magnitude 
of Iab = 28 for point source, and Ia b = 26 for galaxies of 
diameter ~ 0.5" (|Scoville et aLlbOOTJ ). This is much deeper 
than SDSS. The specific COSMOS catalogue used here is 
the Zurich Structure and Morphology catalogue!, whi c h con - 
tain s the measurements presented in IScarlata et al.l (|2007l ) 
and ISargent etahl (120071 ) . 

We further excluded COSMOS and SDSS objects which 
are outside the region common to both the COSMOS and 
the SDSS masks. The SDSS mask used here is made up of the 
"spherical polygons" provided on the NYU-VAGC website, 
and is the same as used in the main body of our paper when 
evaluating and correcting for edge effects. The COSMOS 
mask was kindly provided by Xu Kong (private communi- 
cation). Objects in the Zurich COSMOS catalogue which 
have flag ACS.CLEAN=0, ACS^U.CLASS=3 or JUNK- 
FLAG=1 are also excluded. We consider objects (galaxies 
or stars) to be matched in COSMOS and SDSS if their an- 
gular separation is smaller than 1.5". 

Figure IX2I is a scatter plot comparing the COSMOS and 
SDSS i-band magnitudes of matched objects. For about 4% 
of SDSS objects with r < 21 there are two matched candi- 
dates in COSMOS. This reflects the much higher resolution 
of the COSMOS data which allows substantially better sepa- 
ration of close pairs of images, double stars, merging galaxies 
and chance superpositions, than is possible with SDSS. In 
such cases, we sum the flux of all the matched COSMOS 
objects and compare to that of the single SDSS object. 

The red straight line is a reference with icosmos ~ Udss 
to guide the eye. In general, i-band magnitudes in COS- 
MOS and SDSS agree with each other quite well, with a 
bigger scatter for fainter objects. SDSS r-band magnitudes 
of r = 19, 20, 20.5 and 21 correspond roughly to i-band mag- 
nitudes, i = 18.6, 19.6, 20.1 and 20.6. For these SDSS values 
the scatter in COSMOS magnitudes in figure El is 0.277, 
0.293, 0.360 and 0.367 magnitudes, respectively. 

The histograms at the top and right of figure IA2I show 
the fraction of objects in each survey which is unmatched in 
the other as a function of i-band magnitude. It is expected 



* We use 150.7 > ra > 149.5 and 2.72 > dec > 1.65 to exclude 

the irregular edge of the full survey. 

^ http://irsa.ipac.caltech.edu/data/COSMOS/tablcs/inorphology/ 



that SDSS objects should be lost at the faint end, because 
COSMOS is much deeper. At the bright end, there are a 
few spikes, but the total number of objects at the bright 
end is small. For i < 20.5, we fail to match about 1% of 
SDSS objects in the COSMOS catalogue. This is because 
the COSMOS mask we are using is not identical to that used 
for the Zurich catalogue. In particular, it eliminates smaller 
regions around bright stars so we find a few SDSS objects 
in such regions which are not in the Zurich catalogue. 

At i ~ 20.5, there are also about 1.5% of COSMOS ob- 
jects which are unmatched in SDSS. This is only 71 objects 
in total, and we have looked at both SDSS and COSMOS 
images of these fields. We found that about 31 of these 71 
objects failed to match because the angular separation be- 
tween their centroid coordinates was bigger than 1.5". In 
some cases the objects are quite extended or have a close 
neighbour. Given the relatively poor seeing in SDSS, the 
determination of the coordinates of faint objects can be un- 
certain and so can result in failure to match to COSMOS. 

On the other hand, about 40 of the 71 objects really do 
not have a match in SDSS. Of these, about 1/3 do indeed 
have an object at or close to their coordinates in the COS- 
MOS image which SDSS failed to recognize. The remaining 
2/3 do not, however, have an object at or close to their coor- 
dinates even in the COSMOS image. The Zurich Structure 
and Morphology catalogue al so provides a cross- match with 
the ground based catalogue of lCapak et al.l (|2007l ). Using the 
cross-match information provided, we found 15 o ut of these 
40 objects do have a match in the lCapak et al.l (|2007l ) cat- 
alogue, with 12 having angular separation between the two 
matched objects bigger than 1.5". The remaining obj e cts do 
not have a match in the ground-based ICapak et al.l (120071 1 
catalogue. We are not clear what produced these seemingly 
spurious objects in the COSMOS catalogue, but the fraction 
is indeed small (< 1%). 

The overall conclusion of this exercise is that incom- 
pleteness in the SDSS/DR8 photometric catalogue is at or 
below the 1% level at the apparent magnitude limit relevant 
for our analysis. This is small enough to be negligible. 

A3 Star-Galaxy Separation in SDSS 

SDSS carries out star-galaxy separation using the mea- 
sured difference between psf (Point Spread Function) 
and cmodel (Composite Model) magnitudes. If psfMag- 
cmodelMag > 0.145, an object is classified as a galaxy, oth- 
erwise it is called a star. The quality of SDSS star-galaxy 
separation is thus unavoidably dependent on observational 
conditions such as seeing and sky background, and is also 
expected to be a function of galaxy apparent magnitude and 
angular size. On the other hand, COSMOS should give much 
more reliable star-galaxy separation because of its substan- 
tially greater depth, its much smaller point-spread function 
and its more uniform observing conditions. In this section 
we again use COSMOS data to evaluate the quality of star- 
galaxy separation in the SDSS. 

We take galaxies in the COSMOS field and cross-match 
them to SDSS images, both galaxies and stars. In the match, 
we again used objects with 150.7 > ra > 149.5 and 2.72 > 
dec > 1.65 in order to avoid the COSMOS edge, and we 
again excluded objects which are obscured by the mask of 
either survey. 
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Figure A3. Top left: The fraction of COSMOS galaxies which 
have a matched SDSS galaxy (as opposed to an SDSS "star" ) as 
a function of r-band apparent magnitude in SDSS. : Top right: 
The fraction of COSMOS galaxies which have a matched SDSS 
galaxy with r < 21, reported as a function of observed SDSS g — r 
colour. Bottom left: Similar to the top right panel, but showing 
the matched fraction as a function of angular size as measured in 
COSMOS, and again limited to objects with r < 21 in the SDSS 
catalogue. 



For COSMOS galaxies that have a match in SDSS (ei- 
ther a "galaxy" or a "star"), we show in figure [XS] the frac- 
tion which are classified as galaxies by SDSS as a function 
of SDSS r-band apparent magnitude, SDSS g — r col our and 
COS MOS GIM2D (Galaxy Image Two-Dimensional ISimardI 
1 19981 ) half-light radius. Th e COSMOS G IM2D half-light ra- 
dius is obtained by fitting alSersid (|l968l ) model to the galaxy 
image (|Sargent et al.ll2007h . 

In the top left panel of figure IA3I the fraction of 
"correctly-classified" SDSS galaxies is close to 1 at r < 20 if 
we take COSMOS as reference, and drops significantly be- 
yond r = 21.5. At r = 21, the fraction is about 95%. There 
is some discussion about star-galaxy separation on the offi- 
cial SDSS websitqij. SDSS/DRl was compared against the 
COMBO-17 survey (Classifying Objects by Medium-Band 
Observations in 17 Filters). This much older completeness 
curve also drops significantly beyond r — 21.5, and the frac- 
tion at r = 21 is about 95%, in excellent agreement with our 
COSMOS test here. 

We now consider all COSMOS galaxies with matched 
SDSS apparent magnitude r < 21, and we investigate classi- 
fication success as a function of colour and angular size. The 
top right panel of figure IA3I suggests that blue galaxies are 
more likely to be misclassified than red ones, although the 
correct fraction is always above 90%. Only about 1% of the 
reddest objects are misclassified. We note that the number 
of bluish objects {g — r < 0.7) is only about one-fifth of the 
number of red objects. 

In the bottom panel of figure E3l we report galaxy clas- 
sification success as a function of angular size (Sersic half- 
light radius from COSMOS in units of arcseconds.) The frac- 
tion starts to drop for satellite half-light radii smaller than 
0.4", and is below 90% at 0.2", though never dropping below 
70%. We emphasize that fewer than 100 objects in COSMOS 



■'-'' http://www.sdss3.org/dr8/imaging/othcr_info.php 
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Figure A4. Tests of our background subtraction algorithms us- 
ing stellar mass functions for satellites of isolated primaries in 
the Gil simulation separated into five ranges of primary stellar 
mass as indicated by colour in the legend. Top: Results for iso- 
lated primaries defined in a projection of the full simulation as 
described in the main text and shown previously in figure [6] (the 
solid curves) are compared with the "true" stellar mass functions 
defined as the mean count of satellites within 300 (or 170) kpc 
of these same primaries in 3D (dashed curves). Bottom: Results 
for isolated primaries defined in a full lightconc mock catalogue 
of SDSS/DR8 using exactly the same criteria as for the real SDSS 
surveys and estimated using the identical code used for the SDSS 
data (solid curves). The dashed curves repeat the "true" answer 
from the upper panel. 



have re < 0.2", which is ~ 2% of all galaxies with r < 21. 
Only about 300 have r^ < 0.4". 

In summary, these results together imply that not more 
than about 1% of galaxies brighter than r = 21 are misclas- 
sified as stars by SDSS, which is a remarkably small number 
for a 2.5m telescope working on a less than perfect site. They 
also imply that misclassification is a negligible problem for 
our purposes. 



A4 Background Subtraction Tests 

All our satellite abundance measurements depend critically 
on our ability to construct an unbiased estimate of the ex- 
pected count of unrelated foreground and background galax- 
ies projected close to each primary. This background count 
often exceeds the count of true satellites by a substantial fac- 
tor, and the ratio of the two depends strongly on distance 
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from the primary and on the magnitudes of the primary and 
sateUite galaxies. As a first test of our procedures, we ran- 
domized the positions of aU galaxies in our SDSS/DR8 pho- 
tometric catalogue within the SDSS footprint, while retain- 
ing all their other properties. We then repeated the analysis 
of luminosity and stellar mass functions which led to fig- 
ures [5] and [6] keeping exactly the same sets of spectroscopic 
primary galaxies and using the identical code. This null test 
produced results which were consistent with zero within the 
estimated uncertainties for all satellite luminosity and stel- 
lar mass functions and around primaries in all five stellar 
mass ranges. 

For our next tests, we used the fact that for the Gil 
simulation we know the true properties of the satellite distri- 
bution. We took all central galaxies in the model with stellar 
mass in the five mass ranges we use for our SDSS primaries 
and we compiled luminosity and stellar mass functions for 
their true satellites, defined to be companions within 300 kpc 
of the primary in 3-D. We then compared with the functions 
given in the right panels of figures [5] and [6] which were es- 
timated in projection after correction for background con- 
tamination. The result for the stellar mass function is shown 
in the upper panel of Fig. IA4I Although the normalizations 
differ slightly, as expected because of the different geom- 
etry of the regions within which satellites are counted in 
the two cases, the shapes of the functions agree closely in 
all cases, demonstrating that our "global" correction does 
indeed remove the background properly from our isolated 
galaxy sample. 

Notice that this does not test aspects of our analy- 
sis concerned with converting counts from the magnitude- 
limited SDSS catalogues to the volume-limited statistics 
which we quote, for example, incompleteness corrections 
due to the SDSS mask, effective volume corrections, and 
K-corrections. We therefore took an all-sky "light-cone" 
m ock catalogue construc ted from the MS as described 
in iHenriaues et all (|2012l 1^. We applied our SDSS mask 
and a magnitude limit of r < 21 in order to create a 
mock SDSS/DR8 catalogue. From this we created a mock 
SDSS/DR7 spectroscopic catalogue including a model for in- 
completeness due to fibre collisions and other observational 
problems. We then defined a sample of primaries using ex- 
actly the same criteria as for the real SDSS data (including 
a simple model of our use of photo-z's to deal with spectro- 
scopic incompleteness when checking for isolation) and we 
estimated satellite luminosity and stellar mass functions us- 
ing the same code as for the real data. The results are shown 
in the lower panel of Fig. IA4I Again agreement with the true 
answer is good within the now somewhat larger error bars. 
This agreement provides a full end-to-end test of our code. 

Nevertheless, for real data, there could be additional 
factors which affect the accuracy of background subtraction. 
For example, the surface density of real galaxies on the sky 
could be modulated by extinction, seeing, sky brightness 
or proximity to other objects (see above), and thus affect 
our estimation of background from the whole survey which 
assumes background galaxies to be distributed uniformly. 



^^ This catalogue is available 

|http://www.mpa-garching.nipg.de/millenniuni[ 



at 



Thus it is important to test the quality of our background 
subtraction directly on the real data. 

To do so, we compare measurements of satellite mass 
functions for projected distances between 300kpc and IMpc 
to measurements for projected distances below 300kpc. In 
order to avoid complications induced by our isola tion crite- 
ria wh ich are imposed at 500 kpc and 1.0 Mpc (see lFall et al.l 
(|l976i ) for an example of the strong effects of such criteria) , 
when doing this test we use all SDSS spectroscopic galaxies 
with r< 16.6 without any further selection. From our SDSS 
mock catalogue we estimate the mean surface densities of 
true satellites in these two regions to differ by factors rang- 
ing between 5.34 and 2.62 from the highest to lowest mass 
primaries. Since the background level should be (and is as- 
sumed to be) the same in the two regions and varies strongly 
with satellite mass, any error in background subtraction is 
expected to show up as a distortion in the shape of the satel- 
lite stellar mass function which will be much more marked 
in the outer region. 

Figure Rsl shows results for this test. Solid curves repre- 
sent satellite mass functions measured within Tp = 300kpc, 
while triangles are measurements for IMpc > r-p > 300kpc. 
We renormalize the IMpc > r^ > 300kpc results and re- 
plot them as dashed curves in order to compare their shape 
with the rp < 300kpc results. The agreement between the 
solid and dashed curves is quite good despite the fact that 
the signal-to-background ratio changes by a factor of 3.5 on 
average. This is both encouraging and interesting. Not only 
does it suggest that our background subtraction is working 
well, the close similarity in shape also demonstrates a signif- 
icant regularity of the observed HOD, since the rp <300kpc 
results are dominated by satellites in the same halo as the 
primary and the IMpc > rp > 300kpc results by companions 
in other haloes, at least for the lower mass primaries. 

It is clearly important to carry out a similar test for the 
sample of isolated galaxies which we analyses in this paper. 
For these objects, we need to examine the accuracy of back- 
ground subtraction within 300 kpc, since this is the outer 
radius to which we count companions and lies well inside 
the smallest radius at which we apply an isolation criterion 
(500 kpc). We therefore split the region within 300 kpc into 
two disjoint subregions by radius. We choose the radius sep- 
arating the two subregions so that they would be expected 
to contain equal numbers of satellites if these followed the 
projected NFW profile inferred for the dark matter. We list 
in table lAll the radial boundaries we adopt for each of our 
five different primary mass ranges. 

The top panel of figure IA6I shows satellite stellar mass 
functions as in figure |S] but now split to display results sep- 
arately for the two subregions. Solid and dashed curves give 
results for the inner and outer regions, respectively. Their 
normalizations do not agree exactly, showing that the radial 
profiles are not consistent with the NFW profiles we adopted 
when setting the boundaries. For primaries in the two high- 
est stellar mass ranges (the red and green curves) the ampli- 
tude in the outer annulus is higher than expected relative to 
that near the centre. For primaries in the two lowest stellar 
mass ranges (the cyan and pink curves) the opposite appears 
true, although the measurements are quite noisy. Thus, the 
inner radial distribution of satellites is shallower than NFW 
for massive primaries, and it could be steeper than NFW 
through the main body of the halo for low mass primaries. 
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Table Al. Boundaries of the subregions used for testing baekground subtraetion. 



Range in primary log Af * /Mq 
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Figure A5. Stellar mass functions for companions projected 
within 300kpc (solid curves) and in the range IMpc > 
rp > 300kpc (triangles) of primary galaxies in five disjoint stellar 
mass ranges. The primaries in this plot are spectroscopic galax- 
ies with r < 16.6 but with no isolation criteria applied. Curves 
and symbols of different colour correspond to different stellar 
mass ranges for the primaries, as indicated by the legend. Dashed 
curves show the IMpc > rp > 300kpc results renormalized in am- 
plitude in order to compare their shape to the < 300kpc results. 



W e will leave deta iled discussion of satellite radial profiles 
to I Wang fc White! which is in preparation. 

Despite the difference in amplitude, the dashed and 
solid red, green and blue curves agree in shape quite well, 
except possibly at the lowest masses where statistics are 
relatively poor. The dashed cyan and pink curves are very 
noisy, but within their uncertainties they are also consistent 
in shape with the solid curves. 

The bottom panel of figure IA6I shows the ratio of the 
mean number of background galaxies to the mean number 
of satellites as a function of satellite stellar mass for each 
of the primary samples shown in the top panel. Solid and 
dashed curves again refer to the inner and outer regions. 
For the most massive primaries (the red curves) satellites 
substantially outnumber background galaxies at high mass, 
and are comparable in number at low mass. For the green 
curves, the background count already exceeds the number 
of satellites at almost all masses, while for the blue, cyan 
and pink curves, it is substantially larger, particularly for 
lower mass satellites and in the outer annulus. For these 
lower mass primaries the ratio of background to satellites 
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Figure A6. Top: Stellar mass functions for satellites of isolated 
primaries separated into five ranges of primary stellar mass as 
indicated by colour in the legend. The results shown before in fig- 
ure[6]are split into disjoint inner (solid curves) and outer (dashed 
curves) regions as listed in table EH Bottom: The ratio between 
the mean numbers of background and of satellite galaxies as a 
function of satellite stellar mass. Curves with different colour and 
line type have the same meaning as in the top panel. 



in the outer annulus is almost an order of magnitude higher 
than in the inner region. For the blue curves, where the mass 
functions have high signal-to-noise over the widest satellite 
mass range, the background-to-satellite ratio ranges from 
unity to about 5.5 in the inner region and from 4 to 20 in 
the outer annulus. Give the large values of these corrections 
and their strong variation with satellite mass and radius, it 
is gratifying that the corresponding stellar mass functions 
agree in shape as well as they do. 

A5 Other Observational Factors 

In this section we investigate whether other observational 
factors such as Galactic latitude, extinction, seeing and sky 
brightness introduce systematic effects in our results. We 
start with our sample of isolated primaries with 11.1 > 



© 2011 RAS, MNRAS 000, [T]-?? 



28 Wang et al. 



log M*/M0 > 10.8 and subdivide it into pairs of disjoint 
approximately equal subsamples based on Galactic latitude, 
r-band extinction, seeing, sky brightness at the time of ob- 
servation of each object and redshift. In each case we esti- 
mate satellite stellar mass functions for the two subsamples 
separately and then compare the results. 

In the top two panels of figure IA7I we compare stellar 
mass functions for primaries with Galactic latitude b > 60° 
and with b < 60° as indicated by the legends. (Numbers in 
brackets give the mean Galactic latitudes of the two sub- 
samples.) The two panels differ in that at top left the back- 
ground was estimated globally using the survey as a whole 
and then applied in the same way to primaries in each sub- 
sample, while at top right the background for each subsam- 
ple was estimated separately from its own half of the SDSS 
footprint. A clear difference at low mass is seen between the 
two stellar mass functions in the top left panel, but is al- 
most entirely eliminated in the top right panel. This reflects 
a difference in mean surface density of about 3.5% between 
the two regions for galaxies with (extinction-corrected) ap- 
parent magnitude in the range 21 > r > 20, with the low 
Galactic latitude region having lower mean surface density. 

It is unclear what causes this problem. Clearly there are 
more stars and greater dust optical depths at low Galac- 
tic latitude. However, misclassification of stars as galaxies 
would increase the apparent galaxy density at low latitudes 
and there is no obvious reason why galaxies would be more 
likely to be misclassified as stars at low latitude. Moreover, 
unless the COSMOS field analyzed above is atypical, star- 
galaxy separation in the SDSS is too accurate for errors to 
produce such a change in mean "galaxy" density. Effects 
due to inappropriate correction for reddening also seem to 
be excluded by the explicit extinction test we discuss next. 

The central and lower left panels of figure I A7I show that 
our satellite stellar mass function results do not differ sig- 
nificantly for regions observed through greater extinction, 
in poorer seeing, or against a brighter sky. In all three cases 
the two halves of the sample give results which appear in 
good agreement. Only for the lowest mass objects where te 
uncertainties become substantial is there any hint of a dis- 
agreement between the red and the black curves, and even 
here no consistent pattern emerges. Thus none of the obser- 
vational factors discussed here appears to significantly bias 
our results. 



A6 A bias with redshift? 

In the bottom right panel of figure IXtI we split the primary 
sample according to redshift. There is a significant difference 
in satellite stellar mass function between two resulting sub- 
samples, corresponding to a difference in amplitude by about 
a factor of 1.5 with the nearer primaries having more satel- 
lites. The mean redshifts of the two subsamples are 0.0379 
and 0.0713, and we have checked in the Gil simulation that 
negligible evolution in satellite abundance is expected over 
such a small redshift interval. Further, the distributions of 
stellar mass and intrinsic colour are indistinguishable for 
the two sets of primaries, so a significant "cosmic variance" 
difference between the samples appears excluded. The dis- 
crepancy cannot be a consequence of our isolation criteria, 
since it persists for stellar mass functions like those of fig- 
ure lA6l which count companions around spectroscopic galax- 
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Figure A7. Stellar mass functions for isolated primaries with 
11.1 > logAIir/MQ > 10.8, split into subsamples of approximately 
equal size based on Galactic latitude, r-band extinction, seeing, 
sky brightness and redshift, as indicated by the legends. The top 
two panels show the same subsamples split according to Galactic 
latitude but corrected using a global background estimate on the 
left and a background estimated separately for each subsample 
on the right. Numbers in the parentheses following the legends in 
each panel give mean values of the separating quantity for each 
of the two subsamples. 



ies without imposed isolation criteria. An identical analysis 
of isolated galaxies from the mock SDSS survey we con- 
structed from the MS light-cone finds no difference between 
the satellite stellar mass functions when split by redshift in 
the same way. Thus the discrepancy must reflect a property 
of the real data which is not included in the simulation and 
in construction of the mock catalogue. 

In order to search for clues to the origin of this 
bias, we compare the distributions of rest-frame colour and 
size in these two redshift bins for satellites with 9.5 > 
log M^/Mq > 8.75 which are projected within 300 kpc of 
isolated primaries with 11.1 > logM^/M© > 10.8. Results 
are presented in figure IXSl From the satellite colour distri- 
butions in the left panel, it appears that the higher redshift 
subsample is missing predominantly red satellites. The right 
panel indicates, however, that these "missing" galaxies are 
similar in size to the remaining satellites. The colours, angu- 
lar sizes and surface brightnesses of such objects are in the 
range where SDSS observations detect and correctly classify 
galaxies without difficulty. Hence there is no obvious reason 
why they should be missed if, in reality, they are present. 

Currently, we have no explanation for this apparent de- 
pendence on redshift. We have verified that it appears in 
very similar form for other primary mass ranges. Given that 
the luminosity/stellar mass functions presented in the main 
body of this paper are a combination of counts around pri- 
maries at different redshifts, and that higher redshift pri- 
maries contribute to the counts of brighter/more massive 
satellites, this apparent decrease in satellite abundance with 
redshift will result in measurement of a steeper slope for the 
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Figure A8. Left: Rcst-framc colour distributions, normalized 
per primary, for satellites with 9.5 > log M*/Mq > 8.75 projected 
within 300 kpc of isolated primaries with 11.1 > log M^/A/q > 
10.8 and split according to redshift at 2 = 0.05. Right: Size 
distributions, also normalized per primary, for the same two sets 
of satellites. 



full sample as compared to either the black or the red curve 
in the bottom right panel of figure IA7I The measured low- 
mass slopes for these two subsamples are -1.444 (black) and 
-1.390 (red). The corresponding low-mass slope for the stel- 
lar mass function in figure [6] is -1.455 (see table [l]), which 
is indeed steeper, as expected. The difference is less than 
5%, however, so this residual unidentified systematic should 
have only a small effect on our conclusions. 
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